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ABSTRACT 

The  explosion  limits  associated  with  a  given  sample  of  an  unstable 
compound  are  generally  taken  to  be  defined  by  the  sample’s  temperature, 
pressure,  composition,  and  size.  In  the  presence  of  steady-state  or 
pulsed  radiative  fields,  these  explosion  limits  are  necessarily  modified. 
Further,  photochemically  significant  radiation  (e.g.,  U.V.)  may  play  a 
markedly  different  role  in  this  modification  than  does  purely  thermal 
radiation  (I.R.).  Examination  is  made  of  the  role  of  both  steady-state 
and  pulsed  radiative  fields  in  the  explosive  behavior  of  unstable  com¬ 
pounds.  The  analytic  considerations  aie  illustrated  for  the  case  of  the 
mixed  hydrazine  family  of  reactants  as  well  as  for  ozone. 
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NOMENCLATURE 

c  reactive  intermediate  in  the  decomposition  of  N^H^  and  related 
compounds  (see  Appendix  I) ,  moles /cm^ 

3 

0-atom  concentration,  moles /cm 

1  3 

concentration  of  [  D]  atomic  oxygen,  moles/cm 

a  reactive  intermediate  In  the  decomposition  of  and  related 

compounds  (see  Appendix  I),  moles/cm^ 

3  -  3 

concentration  of  [  diatomic  oxygen,  moles /cm 

concentration  of  a  reactive  intermediate  in  the  ordinary  quasi- 

3 

steady  state,  moles/cm 

concentration  of  a  reactive  Intermediate  in  the  photochemical 

3 

quasi-steady  state,  moles /cm 
heat  capacity,  cal/gm/°K 

characteristic  dimension  (cylinder  diameter)  for  the  sample  of 
unstable  compound,  cm 

diffusion  coefficient  for  species  A  and  the  remainder  of  the 
2/ 

gas,  cm  /s 

diffusion  coefficient  for  species  B  and  the  remainder  of  the 

2, 

gas,  cm  /s 

dissociation  energy,  calories/mole 

energy  of  electronically  excited  oxygen  atoms,  calories/molc 

energy  per  mole  of  photochemical ly  significant  photons,  calories/ 
mole 

3 

reactant  concentration,  tooles/cm 
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3 

■  initial  value  of  reactant  concentration,  moles/ cm 

3 

«  volumetric  radiative  energy  source  function,  calories/cm  /a 

*  Planck's  constant 

■  heat  of  reaction  for  the  ith  kinetic  step,  where  i  ■  1,2,  etc., 
calories/mole 

*«  rate  constant  for  the  jLth  kinetic  step,  where  i  ■  1,  2,  etc. 

■  rate  constant  for  the  photodissociative  step,  s  ^ 

3 

**  heat  loss  rate,  cal/cm  /a 

3 

«  volumetric  (average)  conduction  heat  loss  rate,  calcries/cm  fa 

3 

®  volumetric  radiative  energy  loss  rata,  calories /cm  /a 

3 

*  third  body  concentration,  moles/cm 

*  number  of  moles  of  reactant  decomposed  per  unit  volume  by  a 

3 

short  radiative  pulse,  moles/cu 


[0]  “  concentration  of  ground-state  atomic  oxygen,  moles/cmJ 
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*  concentration  of  electronically  excited  diatomic  oxygen  [  E  ] , 

3  U 
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3 

*  a  product  molecule,  moles/cm' 

3 

“  heat  release  rate,  calories /cm  / a 

■  reaction  rate  for  the  ith  kinetic  step,  where  i  ■  1,  2,  etc., 

3 
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-  reaction  rate  for  the  photodissociative  kinetic  step,  moles/cm  /a 
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T2  ■  unstable  steady-state  temperature  for  the  unstable  medium,  °K 

t  «  critical  temperature  for  explosion  of  the  unstable  medium,  °K 

t  m  environment  (bath)  temperature  for  the  unstable  medium,  °K 
b 

T  »  critical  bath  temperature  for  explosion,  °K 
be 

T  »  temperature  of  the  quasi-steady  sample  of  an  unstable  compound,  K 

s 

X  m  temperature  of  the  sample  of  an  unstable  compound  under  photo— 

PS  *  0 

chemical  quasi-steady  state  conditons,  K 

AB  *  instantaneous  concentration  rise  due  to  an  arbitrarily  short 

3 

radiative  pulse,  moles/cm 

AT  *  instantaneous  temperature  rise  due  to  an  arbitrarily  short 
radiative  pulse,  °K 

X  •  thermal  conductivity,  calories/cm/s/°K 

v  *  frequency,  8  1 

p  »  density,  gm/cm^ 

MMH  «  Hanomethylhydrazine 

UBMH  *  unsymetrical  dimethylhydrazine 
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I.  INTRODUCTION 

&r.  unstable  compound  in  a  given  apparatus  is  characterized  by 
pressure-temperature-composition  regimes  which  correspond  to  explosive 
or  non-explosive  conditions.  These  are  the  well  known  explosion  limits. 
Implicit  in  this  experimental  definition  of  an  explosion  limit  is  the 
absence  of  an  external  radiative  field  which  could  influence  either  the 
thermal  or  the  chemical  properties  of  the  potentially  explosive  system. 
Where  such  a  radiative  field  exists,  its  effects  on  the  explosive  proper¬ 
ties  of  flie  unstable  compound  must  he  considered.  Such  consideration  is 
of  more  than  pure  academic  interest.  For  example,  rocket  engine  environ¬ 
ments  generally  contain  strong  sources  of  radiation  which  can  induce  both 
thermal  and  photochemical  effects  in  an  unreacted  unstable  compound;  such 
as:  hydrazine,  any  member  of  the  mixed  hydrazine  family  of  reactants, 
ozone,  ozone-halogen  reactants,  and  many  others.  The  rocket  engine 
radiative  environment  experienced  by  a  sample  of  an  injected  unstable 
compound  corresponds  to  an  intense  radiative  transient.  Properties  of 
this  radiation  field  (together  with  the  other  engine  environmental 
features)  are  required  to  prescribe  the  explosion  conditions  for  the 
unstable  compound. 

Another  case  of  interest  is  that  associated  with  the  well  known 
flash  photolysis  experiments.  Here  again,  a  transient  radiative  field, 
containing  both  thermally  and  photochemically  significant  radiation,  plays 
a  role  in  the  explosion  phenomenon.  Another  case  of  practical  interest 
is  exemplified  by  a  given  sample  of  an  unstable  compound  subjected  to  a 
far  less  Intense  but  steady  radiative  field.  Such  a  sample  has  its  normal 
explosion  limits  affected  by  the  radiative  field. 

It  is  the  purpose  of  this  study  to  exsmine  the  various  ways  in 
which  photochemically  and  thermally  significant  radiation,  both  transient 
and  steady,  can  affect  explosion  limits  of  unstable  materials. 
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U.  RADIATIVE  EFFECTS  ON  EXPLOSIONS 

Both  energy  and  mass  conservation  equations  nay  be  utilized  to 

prescribe  explosive  instability  features  which  are  unobtainable*  or  less 

2  3 

fully  prescribed  by  classical  consideration  *  of  only  one  equation  or 
the  other. 

For  ultraviolet  absorption  by  an  unstable  sample  (such  as  ozone  or 
hydrazine) ,  the  primary  photochemical  act  involves  creation  of  a  predisso¬ 
ciated  state,  followed  by  decomposition  of  the  absorbing  molecule  and 
partial  thermallzatlon  of  the  absorbed  photon's  energy.  The  thermallzed 
portion  of  this  energy  source  must  be  accounted  for  in  the  energy  conserva¬ 
tion  equation.  The  remainder  of  the  energy  is  used  to  generate  disso¬ 
ciated  molecular  fragments,  and  corresponds  to  a  kinetic  source  term 
that  must  be  included  in  the  mass  conservation  (kinetic)  equations.  Thus, 
the  energy  equation  is  written: 

1  ■  *■=,  £  -  +  I  Bihl  +  <»■ 

Additionally,  the  kinetic  rate  equations  are  coupled  to  equation  (1) 
This  coupling  is  illustrated  below  with  the  simplified  kinetics  as  often 

4 

considered  for  1002  ozone  or  hydrazine  decomposition: 


F  +  h\H-v  B  +  P  (2) 

kl 

F  +  M  B  +  P  +  M  (3) 

k2 

F  +  B  2P  ,  (4) 

where  the  rates  of  the  individual  steps  are  given  by  the  equations 

Rv  -  \(IF]  (5) 
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(6) 


R2  -  k2[F][B]  (7) 

etc.-, 

i 

and  where-  rate*-  for-  the-  individual  species  take  the  form: 

.  o/r  -  *v  -  b1  -  r2  ±  .  (8) 

-  DbV2B  +  Rv  +  Bl  -  R2  +  .  (9) 

etc;. 

When  the  approximation  of  spaticl  homogeneity  is  permitted,  equations 
(1)  ,  (8),.  and  (9)  become:- 

pcp  dk  "  E  ftlhl  +  3th  «°> 

and- 

am .  ..,v .  Bi .  *2  (11) 

f 

”dtf  “  %  +  *1  “  R2. 

■etc;-*,  where. .-all;  values :  are  spatially  averaged. 

^nder.  some,  circumstances  ,  the  explosion:  limits  of  such  a  system 
may ’be  studied  by  the  method  of  the  phase  plane,  or  by  modified  (classical 
methods..  Where  characterization  of  the/chemical  kinetic  scheme  is  more 
complex,. cor  transient  radiative-  fluxes  are;  to  be  considered,  direct 
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integration  of  the  coupled  conservation  equations  may  be  necessary.  The 
analytic  approach  to  be  taken  differs  with  the  nature  and  complexity  of 
the  specific  case  under  consideration. 

It  is  necessary  to  reemphasize  the  fact  that  the  presence  of  an 
ultraviolet  radiative  field  generally  implies  thermal  as  well  as  photo** 
chemical  sources  in  the  unstable  medium.  The  relative  significance  of 
each  type  of  source  depends  on  many  factors,  particularly  the  chemical 
nature  of  the  irradiated  reactant.  In  general,  however,  analysis  of 
photolysis  requires  that  thermal,  photokinetic,  and  the  usual  chemical 
kinetic  effects  all  be  considered. 
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III.  SPECIAL  CASES  INVOLVINS  RADIATIVE  EFFECTS  ON  EXPLOSIONS 

It  is  expected  that  analysis  of  the  explosive  behavior  of  a  photo- 
chemlcally  active  unstable  system  will  depend  upon  such  features  as  the 
optical  opacity  of  the  medium  as  well  as  on  the  spectral  and  temporal 
structure  of  the  incident  radiative  flux  density.  We  here  delineate 
several  limiting  cases  and  indicate  the  kinds  of  analyses  that  may  be 
utilized  in  examining  each  case. 

For  the  purposes  of  this  discussion,  we  assume  the  sample  is  in  a 
bath  of  temperature  and  that,  in  the  absence  of  radiation,  explosion 
does  not  occur. 


CASE  (A); 


This  case  may  be  analyzed  by  the  method  of  the  phase  plane.  If 
is  the  energy  of  an  absorbed  photon  and  (Ey  -  E^)  the  thermalized  portion 
of  this  energy,  then,  the  instantaneous  initial  temperature  rise  of  the 
system  is  given  by: 


For  a  reaction  of  the  form 

F  +  hv  ■+*  2B  +  (Ev  -  Ep) 


(14) 


where  B  is  some  reactive  molecular  fragment  produced  by  photodissociation 
of  F,  there  is  an  "instantaneous"  rise  in  the  value  of  B,  due  to  the 
flash: 

AB  ■  2n,  (moles/cm^).  (15) 
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In  order  to  examine  the  stability  i'T  the  flashed  system  one 
examines  the  new  state  coordinate  (T^  +  at,,  +  AB)  on  the  phase  plane, 
where  (T^,  B^)  arc  the  state  coordinates  prior  to  the  flash.  If 
(T^  +  AT,  +  AB)  falls  on  an  unstable  trajectory,  explosion  follows. 

CASE  (B) :  The  optically  thin  sample  is  subjected  to  constant,  steady- 
state  u.v.  flux. 


This  case  may  be  analyzed  by  one  of  several  methods.  The  steady- 
state  u.v.  flux  changes  the  ordinary  (quasi-steady)  state  to  a  photo¬ 
chemical  steady  state,  or  explosion  occurs.  Thus,  if  (Tc,  B  )  are  the 

b  b 

space  average  values  of  T  and  B,  prior  to  irradiation,  they  become  T^g 
and  shortly  after  irradiation  is  started.  Here  T^g  >  Tg  and  Bpg  >  Bg. 
Either  classical  methods^*^  or  the  phase  plane  method  may  be  utilized 
(as  appropriate)  in  the  determination  of  the  critical  explosion  conditions. 


CASE  (C) :  The  optically  thin  sample  is  subjected  to  a  nonsteady  u.v. 

flux  that  ultimately  becomes  steady. 


If  the  final  steady  state  flux  corresponds  to  an  explosion  condition 
(Case  B) ,  then  explosion  occurs .  If  the  ultimately  achieved  steady  state 
flux  corresponds  to  a  photochemical  quasi-steady  state,  explosion  may  or 
may  not  occur  prior  to  such  time.  In  general,  integration  of  the  time- 
dependent  equations  [e.g.,  equations  (10),  (11),  (12)3  reveals  whether 
explosion  occurs  or  whether  a  photochemical  quasi-steady  state  is 
achieved. 

CASE  (D) :  The  optically  thick  sample  is  subjected  to  an  arbitrarily 
short,  intense  u.v.  pulse. 

These  initial  conditions  leud  to  inhomogeneous  photoinitiation. 

The  time  and  space  integration  of  the  more  general  conservation  .equations 
[e.g.,  equations  (1),  (8),  and  (9)3  may  be  necessary  to  a  determination 
of  these  effects. 


20 


UNCLASSIFIED 


UNCLASSIFIED 


CASE  (E) :  The  optically  thick  sample  i8  subjected  to  a  steady-state 
u.v.  flux. 

These  conditions  lead  to  the  establishment  of  a  spatially  hon- 
homo^eneous ,  photochemical  quasi-steady  state,  or  to  nonhomogeneous 
photoinitiation.  Explosion  is  taken  to  occur  when  there  exists  an  upper 
limit  on  the  u.v.  radiative  flux  density  which  permits  quasi-steady 
solutions  to  the  nonhomogeneous  conservation  equations  [e.g.,  equations 
(1),  (8),  (9)]. 

In  the  analyses  appropriate  to  homogeneous  explosion  one  requires 
knowledge  of  the  chemical  kinetics  of  the  system  involved.  Additionally, 
the  characterization  of  the  photochemical  steady  state  and  the  applicable 
time-dependent  conservation  equations  are  needed.  For  the  mixed  hydrazine 
family  of  fuels,  these  are  considered  in  Appendix  I. 
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IV.  EXPLOSIVE  EFFECTS  OF  A  RADIATIVE  PULSE  ON  A  QUASI-STEADY  SAMPLE 

OF  AN  UNSTABLE  COMPOUND 

Based  on  these  considerations,  we  now  characterize  the  ways  In  which 
a  thermally  and/or  photochemically  significant  radiative  pulse  can 
explosively  modify  the  state  of  a  quasi-steady  sample  of  an  unstable 
compound.  An  illustrative  phase  plane  diagram  is  shown  in  Fig.  1.  At 
time  equals  zero,  and  in  the  absence  of  radiation,  the  system  is  quasi¬ 
steady  and  the  T~B  values  characterizing  the  system*1 2  are  those  of  the 
stable  nodal  point,  (T^,  B^).  The  unstable  saddle  point  conditions  are 
(T^,  B^).  The  separatrix,  S,  defines  the  boundaries  between  explosive 
and  non-explosive  trajectories.  The  explosive  trajectories  must,  in  a 
real  case,  approach  a  stable  singularity  associated  with  thermodynamic 
equilibrium. 

In  this  phase  plane  diagram,  we  have  indicated  various  ways  (both 
real  and  hypothetical)  in  uiich  the  quasi-stable  system  can  be  radiatively 
exploded: 

1.  A  thermal  pul?<.,  auch  as  3>^,  is  represented- by  a  straight 
vertical  line.  „f  the  energy  of  the  pulse  is  sufficient,  the 
change  in  the  originally  3tabla  system  is  characterized  by  a 
sharp  temperature  Increase  and  the  resultant  unstable 
trajectory  leads  to  explosion.  Physically,  such  a  case  can 
only  be  approximated. 

2.  A  thermal  pulse,  such  as  whose  duration  is  comparable  to 
the  response  time  of  the  system,  must  be  represented  by  a 
curved  line.  Although  the  primary  act  of  photoabsorption 
does  not  result  in  dissociation  of  the  unstable  compound  (and 
consequent  production  of  [B]),  the  finite  duration  of  the 
pulse  requires  that  [B] ,  as  well  as  T,  increase  during  the 
lifetime  of  the  pulse.  Such  photoinitiatlon  of  explosion  can 
be  achieved  with  a  powerful  infrared  pulse  and  corresponds 

to  a  physically  realizable  case. 


B] 


(Moles/cc) 
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Figure  1.  Types  of  Reel  and  Hypothetical  Radiative  Pulses  That  Can 
Explosively  Modify  A.  Quasi-Steady  Sample  of  an  Unstable  Compound 
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A  photodissociative  pulse,  such  as  is  represented  by  a 
straight  horizontal  line.  If  the  energy  of  the  pulse  Is 
sufficient,  the  change  In  the  originally  stable  system  is 
characterized  by  a  large  increase  in  the  reactive  intermediate 
concentration,  and  the  resultant  unstable  trajectory  leads  to 
explosion.  Physically,  such  a  case  can  only  be  approximated. 

A  pulse  such  as  3^,  containing  ultraviolet  and  other  wave¬ 
lengths,  is  represented  by  a  straight,  diagonal  line.  If 
the  energy  of  the  pulse  is  sufficient,  the  change  in  the 
originally  stable  system  is  characterized  by  both  a  temper¬ 
ature  increase  and  an  increase  in  the  reactive  intermediate 
concentration  and  the  resultant  unstable  trajectory  leads  to 
explosion. 

Ultraviolet  radiation,  or  a  combination  of  ultraviolet  and 
longe  i/elength  radiation,  leads  to  both'  photodlssoclatlon 
and  thermalization  of  the  incident  radiation.  Physically, 
however,  the  arbitrarily  short  pulse  time  can  only  be 
approximated. 

A  pulse,  such  as  3^  (again  with  duration  comparable  to  the 
response  time),  containing  ultraviolet  and  other  wavelengths, 
is  represented  by  a  curved  line.  The  primary  act  of  photo- 
absorption  (creation  of  a  pre-dissociat ed  state)  leads  both  to 
photodecomposition  and  to  thermalization  of  the  incident 
photon's  energy.  Beth  T  and  [B]  increase  and,  if  the  proper¬ 
ties  of  the  pulse  are  adequate,  the  change  in  the  originally 
stable  system  leads  to  an  unstable  trajectory  and  explosion. 
Such  photoinitiation  of  explosion  can  be  achieved  with  a 
powerful  ultraviolet  pulse,  or  some  mixture  of  ultraviolet- 
containing  wavelengths.  Such  a  case  is  possible  when  an 
unstable  compound  is  suddenly  injected  into  the  radiative 
environment  of  a  rocket  engine.  This  also  is  the  case 
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characterizing  the  well  known  flash  photolysis  experiments. 

Of  course,  photolysis  which  does  not  lead  to  explosion  is 
represented  by  the  stable  trajectories  which  lie  below  the 
separatrix. 

Of  the  above  five  cases,  (2)  and  (5)  represent  systems  which  can 
be  physically  realized.  There  is  a  continuum  of  possible  (AT,  AB)  values 
associated  with  the  possible  pulses  that  can  just  barely  drive  the  originally 
stable  system  onto  an  unstable  trajectory.  Clearly,  the  energy  of  the 
flash  does  not  in  itself  characterize  photoignition,  for  any  given 
unstable  system.  -Spectral  distribution  plays  a  role.  Also,  there  exists 
an  infinitude  of  unstable  trajectories.  Consequently,  it  is  apparent, 
for  the  case  of  flash  photolysis  experiments,  that  the  experimentally 
observed  (T,  B)  trajectories  seen  after  cessation  of  the  flash,  are  not 
necessarily  independent  of  the  initiator.  It  is  also  true,  however, 
that  as  the  unstable  trajectories  converge  to  the  singularity  associated 
with  the  final  products  in  equilibrium,  that  the  (T,  B)  values  along  the 
trajectories  converge,  in  time. 

Implicit  in  the  representation  presented  in  Fig.  1  is  the  approxi¬ 
mation  that  the  initial  reactant  concentration  has  not  been  significantly 
changed  during  the  time  period  of  the  ultraviolet  flash.  This  may  not 
be  a  suitable  approximation  in  all  cases. 

In  the  next  section,  the  calculated  time  behavior  of  temperature 
and  intermediate  species  concentrations  is  presented  for  a  given  case 
of  ultraviolet  irradiation  and  explosion  of  hydrazine. 
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V.  TMg  PHOTOCHEMICAL  STEADY  STATE  AND  EXPLOSIONS 


2-3 

Classically,  thermal  explosion  theory  considers  the  quasi-steady 
state  of  a  sample  of  an  unstable  compound  to  be  at  a  temperature  (Fig.  2a) 
when  the  heat  release  rate  and  heat  loss  rate  of  a  reaction  are  plotted 
versus  temperature.  T2  corresponds  to  an  unstable  point  for  which  the 
heat  Ids 8  and  heat  generation  rates  are  equal.  Explosion  is  taken  to 
occur  when  and  T2  merge,  as  shown  in  Fig.  2b.  In  modifying  such  a 
representation  for  the  presence  of  a  steady  state  photodissociative  flux, 
the  heat  release  rate  function  must  account  for  both  the  thermal  and 
kinetic  effects  of  the  photoabsorption  events.  Here,  either  a  photochemical 
quasi-steady  state  is  achieved,  or  explosion  occurs  (Case  (B)  of  Section 
III).  Again,  explosion  is  taken  to  occur  when  the  heat  loss  rate  and 
heat  generation  rate  functions  are  just  tangent  (as  in  Fig.  2c).  For  the 
case  of  a  cylindrical  sample  of  hydrazine,  the  kinetic  scheme  discussed 
in  Appendix  I  may  be  used  to  obtain  the  following  values  of  [NH2J  and  T 
at  the  photochemical  steady  state  condition: 
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At  the  explosion  condition,  [NH^,  [NH]g,  and  Tg  depend  upon  the  radia¬ 
tive  field.  Of  course,  the  critical  bath  temperatures  are  not  the  same, 
for  any  two  different  radiative  fields,  each  corresponding  to  an 
explosion  condition. 
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On  the  other  hand,  a  more  general  (phase  plane)  analysis  that 
accounts  for  both  kinetic  and  thermal  considerations  may  be  expected  to 
yield  more  reliable  results. 


Equations  (16)  and;  (17)  may  be.  employed,  together  with'  the  consider¬ 
ation  that  the  nodal  point  andth«  saddle  point  of  the  phase  plane  are 
coincident"*  at  the  explosion  limit  to  deduce  the  explosion  limits  of  a 
cylindrical  sample  of.  irradiated  hydrazine: 
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A  similar  expression  say  be  written  for  the  case  of  ozone. 

The  classical  representation  is  not  suited  for  prescribing  the  time- 
dependent  behavior  of.  a  sample  of  an  unstable  compound  during  ultraviolet 
photolysis.  This  is  particularly  true  when  photolysis  occurs  during  a 
time  period  that  is  comparable  to -the  Induction  time  of  the  reaction  or 
when  a  significant  fraction  of  the  initial  reactant  concentration  is 
consumed  during  the  early  portion  of  the ‘photolysis .  For  such  a  situation 
direct  integration  of  the  time-dependent  equations  is  necessary. 

,  This  is  discussed  in  the  next  sections*  first  for  the  mixed 
hydrazine  family  of  reactants,  then  for  ozone. 


CALCULATED  TEMPORAL  BEHAVIOR  OF  T.  NH.  HH..  N^,  AND  FOR  THE 
ULTRAVIOLET  PHOTOLYSIS  OF  HYDRAZINE 

Employing  the  kinetic  scheme  given  in  Appendix  1  we  examine  worn* 
trajectories  resulting  from  the  ultraviolet  flash  irradiation  of  initially 
quasi-; steady  samples  of  hydrazine.  Similar  calculations  can  be  carried 
out  for  other  members  and  combinations  of  members  of  the  mixed  hydrazine 
family  of  fuels.  The  conservation  equation*-  are; 
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A.  STEADY  STATE  IRRADIATION 

Figure  3  shows  T,  [N^j ,  ~  Uyy,  [Nil],  [NR£)y  andjN^I  as  a 

function  of  tine  for  the  ase  of  photolysis  at  ^250  A,  where  Tb»400°K; 

an  initial  [bLK^]  value  corresponding  to  100%  gaseous  hydrazine  at  a  density 

of  5.69  x  10  •  gm/cm3;  d  «  0.60  cm;  57,100  cal/a»le;  and  a  steady 

state  kv  of  303  a  \  corresponding  to  a  local  ultraviolet  flux  deposition 

rat2  of  6.36  x  103  cal/cm3/s.  The  k  value  goes  from  zero  at  90  micro- 

^  -1 

seconds  to  its  steady  state  value  of  303  s  in  20  microseconds.  Focus¬ 
sing  our  attention  on  the  behavior  of  T  and  [NH^  we  see  that  the  initiating 
process  is  very  much  like  that  indicated  by  3^,  of  Fig.  1.  Cace  an 
explosive  trajectory  Is  reached,  the  behavior  in  the  [T]  vs  [NH  ]  plane 
of  the  system  is  virtually  independent  of  the  (steady  state)  existence 
of  the  ultraviolet  source.  This  has  been  confirmed  by  computations  lit 
which  the  flash  is  cut  off  at  times  such  as  260  microseconds,  270  micro¬ 
seconds,  280  microseconds,  etc.  Ac  substantially  Bhorter  times,  of 
course,  cut-off  of  the  flash  leads  to  a  non-explosive  condition  of  the 
reacting  medium.  Of  particular  interest  is  the  approximately  steady-state 
behavior  of  both  (d/dt)  [N^J  and  [NH2]  in  the  neighborhood  of  200-250 
microseconds,  immediately  before  all  the  derivatives  "blow  up,”  Also  to 
be  noted  is  that  during  the  period  from  100  to  250  microseconds  the 
reaction  process  can  be  regarded  as  quasi-steady  as  far  as  [NH2]  and 
w  are  concerned  and  a  small  but  significant  fraction  of  the 

reactant  is  used  up  before  this  quasi-steady  state  is  destroyed.  Thus, 
the  explosive  behavior  (in  the  neighborhood  of  250-300  microseconds) 
follows  the  destruction  of  the  photochemical  quasi-steady  atate  of  the 
reaction  rate  and  [NH23 .  This  occurs  ±5  the  neighborhood  of  900°K  (see. 

Fig,  3) ,  some  350°K  higher  than  classical  analysis  predicts  (as  discussed 
ia  Section  VI).  The  result?  of  classical  analysis  of  the  photochemical 
steady  state  are  indicated  in  Appendix  ill.  > 

Br  EFFECT  OF  PULSE  DURATION  >" 

In  the  previous  section,  we  discussed  the  irradiation  of  a  N^ 
sample  by  a  square  wave  pulse.  This  feature  is  reflected  by  the,  fact 
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that  *  303  s  ^  throughout  the  time  interval  of  primary  interest. 
Inasmuch  ae  k^  is  proportional  to  the  locct  rate  cf  absorption  of  u.v. 
photons,  a  pulsed  (finite)  radiative  source  corresponds  to  a  k^  that 
goes  from  zero  to  a  finite  value  and  then  back  to  zero  all  within  the  time 
interval  of  interest.  The  effect  of  such  pulses  has  been  investigated. 

In  the  tabulation  given  below,  the  half-life  of  the  sample  is  com¬ 
pared  with  the  time  at  which  the  value  *  303  8  ^  is  suddenly  taken  to 
zero  (pulse  cut-off) . 


Pulse  Cut-off 


Sample  Half-Life 


Microseconds 
Steady  pulse 
294.0 
290.0 
270.0 
240.0 


Microseconds 

295.6 

795.6 

296.4 

335.5 
r  500 


As  the  pulse  is  cut  off  at  shorter  and  shorter  times,  the  sample 
half-life  increases  by  orders  of  magnitude.  Of  course  still  shorter 
pulses  of  the  same  energy  lead  to  half-lives  that  correspond  to  the  "slow 
decomposition"  process.  In  steady-state  analyses,  this  is  the  "non-explosive" 
regime. 

C.  EFFECT  OF  SPECTRAL  PROPERTIES  OF  THE  PULSE 

Given  a  pulse  of  fixed  power  and  duration,  it  is  expected  that 
the  spectral  distribution  of  the  pulse  energy  will  influence  the  temper¬ 
ature-composition  history  of  the  sample.  This  temporal  behavior  depends 
very  greatly  on  the  chemical  nature  of  the  compound  under  consideration. 

For  the  case  of  hydrazine  we  have  performed  calculations  of  these 
trajectories  for  various  spectral  distributions  of  continuous  pulses  at 
constant  total  power  corresponding  to  ky  *  303  s  \ 
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Cases 

Fraction  of 
Pulse  Power 
Thermalized 

Fraction  of 

Pulse  Power  into 
Dissociation 

Sample  N_H, 
Half-Life  4 
(microseconds) 

1  (hypothetical) 

0.00 

1.00 

270 

2  (the  treated  case) 

0.55 

0.45 

295 

3  (a  possible  case) 

0.75 

0.25 

337 

4  (a  possible  case) 

0.90 

0.10 

357 

5  (a  possible  case) 

1.00 

0.00 

373 

For  smaller  values  c?  radiative  flux  density,  the  calculated  spread 
of  sample  half-life  increases  sharply,  reflecting  the  differences  in 
trajectory  times  close  to  the  separatrix  (Fig.  i).  Also  to  be  noted  is 
the  fact  that  Case  5  corresponds  to  a  purely  thermal  pulse  (3^  of  Fig.  1) 
but  that  Cases  1-4  all  lie  in  the  quadrant  associated  with  3?,.  (Fig.  1). 

D.  EFFECT  OF  SAMPLE  DIAMETER 

In  steady-state  explosion  theory,  the  critical  conditions  for 
explosion  depend  on  the  thermal  losses  for  the  system,  and  therefore  on 
the  tube  diameter.  By  the  same  token  (Eq.  (16))  the  rate  of  temperature 
rise  and  the  sample  half-life  of  an  Irradiated  sample  also  depend  on  the 
sample's  diameter.  The  following  results  are  obtained  for  a  constantly 
irradiated  sample  under  conditions  previously  employed: 

Sample  Diameter  Sample  Half-Life 

on  Microseconds 

0.600  295.6 

0.060  296.5 

0.006  485.0 

The  sample  diameter  has  very  little  effect  on  the  sample  half-life 
for  si;.  greater  than  a  few  hundred  microns.  For  explosive  trajectories 
involving  longer  half-lives  (and  smaller  values) ,  the  effect  of  sample 
size  becomes  progressively  greater. 
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VII.  MONOMETHYLHYDRAZINE  AND  ASYMMETRICAL  DIMETHYLHYDRAZINE 

16 


Both  UDMH  and  MMH  are  easier  to  dissociate  than  N.H. , 

2  4 


Consequently 

both  the  heats  of  reaction  and  the  rate  constants  given  for  the  scheme 
(Appendix  I)  must  be  modified.  The  primary  modifications  are  indicated 
below : 


N2H4 

MMH 

UDMH 

*D 

57,100  cal /mole 

51,900  cal/mole 

49,600  cal/mole 

E1 

57,100  cal/mole 

51,900  cal/mole 

49,600  cal/mole 

*5 

18,000  cal/mole 

12,800  cal/mole 

10,500  cal/mole 

<vv 

69,400  cal/mole 

74,600  cal/mole 

76,900  cal/mole 

kl 

.  .19  T  57.1001 

10 

..19  T  51.9001 

i°  «p  [-  -fe-j 

. .19  T  49.6001 

10  “*r  sri 

k5 

io12  «,[-  ugu| 

,„12  r  12.8001 

!0  exp  j-  -Jj-J 

Although  the  optical  extinction  coefficients  for  the  above  three 
reactants  are  approximately  equal,  the  lower  E^  is  the  greater  the  fraction 
of  an  absorbed  u.v.  photon's  energy  which  goes  directly  into  raising  the 
sample's  temperature.  Additionally,  the  lower  is  the  higher  k^  and 
are.  Thus,  from  both  a  thermal  and  photochemical  point  of  view,  UDMH  is 
the  most  reactive  of  the  three  reactants  and  the  least  reactive. 

No  consideration  is  given  here  to  possible  heterogeneous  reactions. 


A  test  of  the  kinetic  scheme  and  rate  constants  was  made  for  the 

case  of  UDMH.  Figure  4  shows  a  comparison  of  our  calculated  explosion 

14 

limits  with  those  observed  by  P.  Gray  and  M.  Spencer.  Agreement 
appears  good,  particularly  when  one  considers  the  unsuitable  geometry 
(short  cylinder)  used  by  Gray  and  Spencer. 
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For  ambient  conditions  similar  to  those  chosen  for  ^H^,  .steady- 
state  explosion  limits  were  calculated  for  UDMH  and  MMH.  For  p  ■  46  atm, 
MMH  (sample  size  of  0.6  cm  diameter)  becomes  thermally  unstable  at  an 
average  sample  temperature  411°K;  for  UDMH,  this  occurs  at  361°K.  As 
noted  earlier,  thermal  instability  does  not  necessarily  assure  a  sharp 
thermal  and  pressure  spike. 

For  the  case  of  MMH  at  a  bath  temperature  of  300°K,  a  steady-state 

-3  2 

local  u.v.  radiative  flux  density  of  4.97  x  10  watt/cm  causes  thermal 

-3  2 

instability.  For  the  case  of  UDMH,  2.41  x  10  watt/cm  causes  thermal 
instability. 

Time-dependent  calculations  show  that  UDMH,  in  a  strong  u.v. 
radiative  field,  gives  a  thermal  spike  in  a  shorter  time  than  MMH  which, 
in  turn,  "explodes"  in  a  shorter  time  than  N^H^.  These  conclusions  are 
implied  by  the  homogeneous  gas  phase  kinetics  noted  previously  and  are 
demonstrated  by  the  calculated  results.  These  are  indicated  below. 

'  _4  3 

For  equal  initial  molar  densities  (1.77  x  10  moles/cm  ) ,  T^  *  40( 
d  ■  0.6  cm,  and  a  k^  ■  303  s  ^  for  t  110  ys,  the  following  results  are 
obtained: 

Sample  Half-Life 
(microseconds) 
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vm.  SOME  SUMMARY  COMMENTS  ON  THE  BEHAVIOR  OF  UNSTABLE  MIXED  HYDRAZINE 

REACTANTS  IN  AN  ENGINE  ENVIRONMENT 

Based  on  our  experimental  observations  as  well  as  our  detailed 
calculations  o£  sample  half-lives,  explosive  trajectories  and  explosion 
limits  for  N^H^,  MMH,  and  UDMH  In  the  presence  and  absence  of  radiative 
fields  (both  U.V.  and  I.R.)  we  make  several  observations: 

1.  Xn  the  absence  of  a  radiative  field  the  explosion  limits  of 

MMH,  UDMH  and  50%-50%  mixtures  of  NgH^-UDMH  occur  at  temperatures 

(<  600°K)  that  are  easily  achieved  in  the  "low  temperature"  regimes  of 
engines.  However,  the  rate  of  temperature  (and  pressure)  rise  in  the 
neighborhood  of  the  explosion  conditions  determines  the  apparent  stability 
of  a  compound  in  a  piece  of  hardware.  Our  experiments  and  calculations 
show  that  this  rate  of  rise,  for  these  compounds,  is  relatively  slow  in 
the  neighborhood  of  600°K.  The  behavior  for  the  case  of  is  typical. 

Just  below  600°K  (Appendix  III)  a  sample  of  ^H^  may  become  explosively 
unstable  (the  steady-state  explosion  condition  is  satisfied).  However, 
a  really  short  sample  half-life  (calculated  to  be  about  50  us)  corresponds 
to  a  gaseous  sample  at  900 °K.  Of  course,  a  reactant  half-life  that  is 
"short"  corresponds  to  a  "short  induct# on  period";  a  'long"  half-life 
corresponds  to  a  "long  induction  period."  Thus,  in  the  absence  of 
additional  external  heat  inputs,  an  explosively  unstable  sample  at  600°K 
(above  the  steady-state  explosion  limits)  has  a  sample  half-life  of  many 
milliseconds  (and  perhaps  much  longer).  As  a  result  of  this  behavior 
(confirmed  both  experimentally  and  analytically) ,  the  engine  designer 
can  usefully  employ  hydrazine  at  temperatures  that  are  higher  than  those 
associated  with  the  explosion  limits. 

2.  In  the  presence  of  a  weak,  steady-state  radiative  field,  the 
explosion  conditions  (associated  with  the  photochemical  steady  state)  are 
achieved  at  a  lower  temperature  than  normal.  Again,  because  of  the  long 
reactant  half-lives  at  the  explosion  condition,  a  reactant  sample  disappears 
slowly  (over  a  period  of  many  milliseconds ,  or  even  seconds  and  minutes) . 
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Rapid  decomposition  does  not  occur  until  the  neighborhood  of  900°K  is 
achieved.  Thus,  the  weak  radiative  field  does  not  really  affect  the 
salient  feature  of  the  system  -  namely  that  the  engine  design  can  usefully 
permit  the  flow  of  the  reactant  in  the  combustion  chamber  at  a  temperature 
above  the  explosion  limit,  for  several  milliseconds,  before  stable  com¬ 
bustion  processes  consume  the  reactant. 

3.  In  the  presence  of  a  strong,  steady-state  radiative  field, 
the  explosion  conditions  (associated  with  the  photochemical  steady  state) 
are  achieved  at  a  still  lower  temperature.  Again,  because  of  the  relatively 
slow  decomposition  behavior  of  at  these  temperatures  (<  600°K) ,  the 
half  life  of  the  sample  can  be  short  (less  than  a  millisecond)  only  if 

the  radiative  field  is  sufficiently  strong  to  drive  the  sample  temperature 
to  some  900°K  in  a  matter  of  a  few  hundred  microseconds.  For  strong 
radiative  fields  (Section  VI),  this  is  indeed  possible.  In  fact,  the 
infrared  radiation  of  an  engine  environment  is  almost  as  effective  as  the 
photochemically  significant  ultraviolet  radiation,  for  the  casd  of  hydrazine 
(Section  VI-C).  (This  is  not  true  for  the  case  of  ozone.)  Inasmuch  as 
rocket  engines  are  strong  sources  of  infrared,  it  is  clear  that  this 
source  of  energy  must  be  included  in  the  engine  designer's  calculation 
of  the  temperature  history  of  a  sample  of  an  unstable  compound. 

4.  It  is  incorrect  for  an  engine  designer  to  estimate  the  temper¬ 
ature  rise  of  a  reactant  sample  (due  to  external  heat  sources,  radiative 
or  other)  without  including  the  self-heating  effects  due  to  the  exothermic 
decomposition  of  the  sample.  Thus,  the  considerations  employed  here,  in 
deducing  explosive  trajectories  and  sample  half-lives,  are  necessary  ones. 
Where  other  than  radiative  terms  play  a  role  in  this  temperature  rise, 
these  new  terms  must  be  included  also  in  the  time  dependent  energy  equation. 

A  particularly  valuable  lesson  to  derive  from  the  behavior  of  this 
class  of  reactants  is  that  the  calculation  of  a  simple  steady-state 
explosion  limit,  in  itself,  is  not  adequate  to  describe  the  stability  and 
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behavior  of  a  reactant  sample  In  an  engine  environment.  Thus,  for  the 
case  of  the  hydrazine  family  of  reactants,  the  "thermal  instability"  of 
the  liquid  does  not  lead  to  a  sharp  pressure  pulse  such  as  that  which  is 
associated  with  a  sample  half-life  of  a  few  hundred  microseconds.  Rather 
this  "thermal  instability"  results  in  an  increased  rate  of  vaporization 
of  the  unstable  liquid  (drop  or  film),  due  to  the  Internal  heat  generation. 
This  occurs  in  a  temperature  regime  where  the  gas  phase  is  not  explosively 
unstable.  Consequently  the  accelerated  vaporized  i  process  (which  may 
give  rise  to  a  "popping"  effect)  acts  as  a  heat  ..ink,  fixing  an  upper 
limit  on  the  temperature  of  the  liquid  sample.  Any  destructively  large 
pressure  pulses  that  may  then  result  must  be  associated  with  the  later 
temperature  history  of  the  resulting  gas  phase  sample.  The  engine 
designer  may  then  determine,  from  the  time  dependent  calculations  (such 
as  those  performed  here),  how  much  time  is  available  for  the  resultant 
gaseous  sample  to  be  mixed  and  burned,  if  the  explosive  behavior  of  the 
reactive  gas  is  to  be  avoided.  For  the  case  of  the  hydrazine  class  of 
reactants,  this  is  easily  achieved  in  the  neighborhood  of  600°K  (sample 
half-lives  of  many  milliseconds)  and  difficult  to  achieve  in  the  neighbor¬ 
hood  of  900°K  (sample  half-lives  of  a  few  hundred  microseconds  or  less) . 
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IX.  CALCULATED  BEHAVIOR  OF  T.  0.  0*,  0 £,  0*,  AND  03  FOR  THE  ULTRAVIOLET 
PHOTOLYSIS  OF  OZONE 

Following  the  analytic  approach  previously  discussed  for  the  time- 
dependent  explosive  behavior  of  the  mixed  hydrazine  family  of  reactants, 
we  consider  the  appropriate  equations  for  the  case  of  ozone.  Employing 
the  kinetic  scheme  discussed  in  Appendix  II  and  setting 

I°3!  -  »•  ‘VlnltLl  ■  V  101  ’  Al’  10,1  ’  V  t°*1  ■  "• 

we  may  write  the  rate  equations  (appropriate  to  high  initial  ozone  con¬ 
centrations)  : 

-  (4r)  ■  k  ?  +  k,  F  F  +  k«A1F  +  LFA  +  k,  F  B  -  f  k,  A  (F  )  (F  -F)  (21) 

\dt  /  v  1  o  21  3  e  4  e26eoo 


Bt  -  ki  F  Fo  -  k2  F  Ai 


“l-  kF  +  k,  F  B  -  k,  F  A  -  k,  A  F 
t/v  4  e  3  e  6eo 


[f<VF>] 


-H-  k„A  F  -  k,B  F  -  krB  F 
t  /  3  e  4  e  5  e  o 


The  energy  equation  yields 


K  (E  -En-E  )  +  L?  h,  +  k,.A.h„ 
v  v  D  e  lol  212 


“  -  —  +  k,A  h«  +  k,  B 

dt  pc  3  e  3  4 

P 


/F\  fl6X(T-Tb) 

A*V.^*»  -  — - 

L  p  J 


m 


F  -Ff A  h, 

O  /  8  6 


Where  assumption  of  a  photochemical  steady  state  is  acceptable,  the 
energy  equation  is: 


pcp!L£ 


as 


'  [('  «)„  ‘  4Bc  +  VWt)]  -  j~ fi'-j  -  R  - 


*  (26) 


and  the  steady  state  concentration  and  decomposition 

[T  F 

*J  -ir 

Tss  2 


rate  equations  are: 

(27) 


[AJ  -v 

•ss 


WoF  +  Wofv^Wol)MJ 


(28) 


(29) 


(30) 


Equations  (26)  to  (30  are  required  for  calculation  of  the  explosive 
stability  of  the  photochemical  steady  state. 


Calculation  of  the  tine-dependent  explosive  behavior  of  a  given 
ozone-containing  sample  requires  consideration  of  Eqs.  (21)  -  (25) 
inclusive. 


Employing  the  kinetics  end  quantum  yield  data  referred  to  in 
Appendix  II,  we  obtain  a  selected  list  of  enthalpies  snd  rates: 


if.%* 


(a)  Enthalpies ;  Ed  -  25,051;  -  101,000;  E^,-  46,000;  hx  *  -  25,051 

h2  -  93,052;  h3  -  2,852;  h4  -  70,949;  h5  *  142,000;  hg  -.71,051  cal /mole. 

(b)  Rates;  k,  -  5  x  1015  exp  [-  24,000/RT];  k,  «  3  x  1013 

exp[-  6000/RT] ;  k^  -  10  exp[-  an  /moie-s;  and  ~  k^  and,  at 

300°K,  k5~  109'0  and  kfi~  105*9. 


A.  PHOTOCHEMICAL  STEADY-STATE  AND  EXPLOSIONS  OF  OZONE 

In  considering  the  explosive  disappearance  of  the  photochemical 
steady  state  of  ozone,  we  note  that  the  temperature  characterizing  the 
rates  of  Eqs.  (27)  -  (30)  is  not  the  bath  temperature,  T^.  Rather,  it 
is  the  temperatures  prescribed  by  simultaneous  solutions  of  (27)  -  (29)  in 
conjunction  with  the  steady  state  assumption  for  the  energy  conservation 
equation: 


• • » 1  dT  . 

q  “  pcp  dt  "  °»  or 


kv[03](Ev  -  ED  -  Ee)  +  kx  t03][M]  hx  +  k2  [03f(0]  h2  +  k3  [03J[C  ]h3j 
+  k4  [03] [0*]h4  +  k5  [0*] [03]h5  +  k6  [03] [M] [0*]h6 


F  (31) 


For  values  of  ky  and  Tfe  low  enough  to  maintain  a  photochemical  steady 
state  that  is  far  from  explosion  conditions,  the  approximation  T^  ~  ^ 
may  be  a  good  one.  In  general,  and  particularly  in  the  neighborhood  of 
an  explosion  radition,  T^  >  T^  nd  the  .approximation  T^  -  T^  is  at  its 
worst. 


42 


UNCLASSIFIED 


Employing  this  thermokinetic  consideration  of  the  photochenical 

if  ^ 

steady  state,  the  steady-state  values  of  [0],  10  ],  [O.^J,  and  are 
deduced  for  a  given  system  (size  and  pressure)  for  a  range  of  bath  temp¬ 
eratures  and  ultraviolet  flux  absorption  rates.  For  a  2.50  cm  diameter 
cylinder  containing  100%  ozone  at  a  pressure  of  0.05  atm,  the  results 
shown  in  Figs.  5-7  are  obtained. 

Examination  of  Figs.  5  and  6  shows  that: 

(a)  (T^  -  Tg)  can  be  a  kinetically  significant  quantity  in  the 
Interpretation  of  a  quantum  yield  experiment  or  an  explosion  limit 
experiment. 

(b)  The  radiative  field  causes  a  decrease  in  the  critical  temper¬ 
atures  for  explosion,  2  and  Tfe  c  from  those  characterizing  the  case 
where  ky  •  0 

Figures  (7a)  to  (7e)  suggest  the  important  roles  played  by  [0] , 

[0  ],  and  [0_]  in  definition  of  the  photochemical  steady  state.  At  lot? 

L  *i 

temperatures  and  .low  k^,  the  heat  release  rate  is  dominated  by  [O^J .  Both 
k^  and  h^  are  large. 

if 

In  the  neighborhood  of  the  explosion  limit,  d[0]/d  k  >>  d[0  ]/dk 

ft  V  V 

>  d(02]/dky.  This  implies  that  simple  two-dimensional  phase-plane 

representation  of  explosion  limits  may  be  adequate  for  representation  of 

the  photochemically  influenced  explosion  limits  of  ozone.  Thus,  for 

constant  k  ,  only  T  and  [0]  need  be  considered  in  the  phase  diagram.  Cf 
^  *  * 

course,  the  roles  of  [0  ]  ana  [0^3  in  determining  the  steady-state  con¬ 
centration  of  [0]  must  be  included. 

B.  EXPLOSIVE  TRAJECTORIES  OF  OZONE 

The  explosive  trajectory  for  a  given  case  of  ozone  irradiation  can 
be  calculated  by  use  of  Eqs.  (21)  -  (25).  This  type  of  calculation  is 
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Figure  6.  Critical  Conditions  for  the  Photochemical  Steady 

State  of  Ozone 
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Figure  7.  Cont'd. 
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shown  for  the  case  of  gaseous  ozone  at  a  pressure  of  1  atmosphere  and  an 

initial  temperature  of  300°K.  A  local  radiative  flux  density  of  100  watts / 
2 

cm  (u.v.)  is  assumed  for  the  6  mm  diameter  cylindrical  sample. 

The  calculated  trajectories  are  shown  in  Figs.  8a  to  8e,  where  the 

'ft  ft 

time-dependent  behavior  of  T,  0^,  0,  0  ,  and  02  are  shown. 

C.  EFFECT  OF  FLUORINE  ADDITIVES  TO  OZONE 

It  may  be  expected  that  the  addition  of  quantities  of  fluorine  to 
ozone  will  introduce  new  chain  carriers,  thereby  sensitizing  the  decom¬ 
position  of  0^.  Thus  reactions  involving  F,  F^ ,  OF 2 ,  etc.  play  a  key 
role  in  reduction  of  the  temperatures  associated  with  ordinary  explosion 
limits  for  ozone.  In  addition,  the  photodissociation  of  F2  and  OF2  (as 
well  as  0^)  in  a  rocket  engine's  radiative  field  suggests  that  trace 
amounts  of  F2  or  OF,,  will  serve  to  sensitize  C>3,  in  a  rocket  engine 
environment. 

On  the  other  hand,  large  concentrations  of  F2  in  ozone  introduce 
a  stable  sensitizer  (F2)  at  the  partial  elimination  of  the  unstable 
reactant  (0^).  Thus,  at  high  concentrations  of  F2  *n  O3  (e. g . ,  where 
[F2J/[03]  >>  1)  the  temperatures  associated  with  the  ordinary  explosion 
limits  for  ozone  must  increase.  Thus,  a  prime  goal  of  a  study  of  explosion 
limits  for  03~F 2  or  03~OF2  fixtures  is  to  define  these  explosion  limits 
for  a  range  of  concentrations  and  radiative  fields.  It  is  clear,  however, 
that  for  high  fluorine  concentrations  (with  ozone),  adequate  explosive 
stability  of  tne  03“F2  sy8tsm  can  be  achieved.  The  detailed  definition 
of  these  explosion  limits  is  a  matter  for  future  study. 
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PURE  OZONE,  G>  Tg  -  300°K,  P  -  1  atm 


(c)  [0]  vs  time 
(Fig.  8  cont’d) 
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X.  SOME  CONCLUDING  REMARKS 

The  roles  of  both  steady  and  non-steady  radiative  fields  in  the 
explosive  behavior  of  unstable  compounds  have  been  examined.  The  relations 
deduced  have  been  employed  in  studying  the  explosive  behavior  of  the 
mixed  hydrazine  family  of  reactants  as  well  as  the  ozone  system. 

On  the  strength  of  these  studies,  it  is  concluded  that: 

(a)  both  hydrazine  and  ozone  systems  have  their  explosion  limit 
condtions  modified  by  the  presence  of  radiative  fields; 

(b)  on  a  per  calorie  basis,  infrared  radiation  is  not  as  effective 

-in  modification  of  the  explosion  limit  as  is  photochemically  significant 

radiation.  This  is  particularly  true  for  the  case  of  ozone,  where  the 

explosion  limit  temperature,  T.  _,  is  greatly  lowered  by  the  photochemically 

£  ^  ft 

induced  concentrations  of  0  and  0^1 

ic)  ozone  is  (explosively)  more  sensitive  to  radiation  the",  is 
hydrazine  5 

(d)  explosive  stability  in  a  rocket  engir-°.  environment  is  of  prime 
interest.  Reactant  half-lives  (ignition  delays)  of  less  than  a  few  milli¬ 
seconds  are  of  major  concern  to  the  engine  designer.  We  have  shown  that 
the  half-life  of  an  explosively  unstable  compound  is  a  very  sensitive 
function  of  the  engine's  radiative  field,  as  well  as  a  function  of  che 
usual  kinetic  and  thermophysical  constants ; 

(e)  the  full  definition  of  the  range  of  conditions  over  which  a 
given  unstable  compound  can  be  employed  in  a  radiative  environment  can 
now  be  made  by  employing  the  deduced  formulations  to  the  deduction  of 
reactant  half-lives  and  explosion  limits,  for  a  range  of  radiative  fields, 
pressures,  temperatures,  and  sample  sizes.  Such  detailed  studies  are 
necessary  before  rational  use  of  ozone  and  ozone-fluorine  propellant 
systems  can  be  made  by  engine  designers. 
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XI.  FUTURE  STUDIES 

Future  studies  of  the  explosive  behavior  of  the  ozone  system  are 
planned  to  include: 

(a)  pressure  as  well  as  temperature  and  composition  trajectories 

(b)  the  effects  of  spectrally  more  complicated  radiative  fields 
(simultaneous  u.v.,  visible,  and  infrared  radiation);  and 

(c)  exploraf!  of  explosive  trajectories  associated  with  longer 
ozone  half-j..  ves  (milliseconds) . 
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APPENDIX  I 

KINETICS  OF  DECOMPOSITION  OF  HYDRAZINE 
AND  RELATED  COMPOUNDS 

6—38 

There  exists  a  substantial  literature  concerning  reaction  rates 
and  photochemical  processes  for  Hydrazine,  Monomethylhydrazine  and  Unsym- 
metrical  Dimethylhydrazine.  For  purposes  of  this  investigation,  we  require 
a  kinetic  scheme  with  knowledge  of  the  rate  constants,  energetics  of  the 
individual  steps,  thermophysical  properties  of  the  reactants  and  products, 
ultraviolet  quantum  yields  and  the  associated  photochemical  extinction 
coefficients.  Although  the  literature  is  not  fully  adequate  to  our 
requirements,  sufficient  data  are  available  to  permit  initial  treatment 
of  the  problem.  Accordingly,  for  the  case  of  hydrazine,  the  scheme  pro¬ 
posed  by  Eberstein  and  Glassman^  is  simplified.  Incorporating  photo- 
dissociative  processes  in  this  modified  scheme  we  write: 

k 

N2H4  +  hv  2NH?  +  (Ev  -  Ed) 

kl 

N0H.  +  M  2NH.  +  M  +  h. 

2  4  2  1 

k2 

N2H4  +  NH2  ->  NH3  +  N2H3  +  h2 

k5 

N2H3  +  M  -*■  NH  +  NH2  +  M  +  h5 

N2Ha  +  NH  NH2  +  N2H3  +  h6 

k7 

N2H3  +  NH2  NH3  +  N2  +  H2  +  h? 

Further  simplifications  may  be  achieved  if  k^  is  taken  to  be  small  (step 

5  is  about  as  endothermic  as  step  1)  and  the  other  k  -  values  modified  to 

reflect  thi«  change.  With  this  modification,  appropriate  schemes  for  MMH 

and  UDMH  decomposition  can  be  represented  in  a  manner  that  parallels  that 

7 

for  NjH^.  Thus,  the  scheme  proposed  by  Cordes  for  UDMH  is  modified: 


t 
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H 

H 


~  (0H3)2  " 


+  NH. 


NNH2  +  M  — ' V  (i 


s)2» 


CH,J  N  +  NH_  +  M 


H 

H 


NNH2  + 


Ha  *  ~H 


HN  +  OH3CH2N  -  NH2 


2b 


NNH2  +  NH2  +  CH0CH„N  -  NH„ 


3  2 


ch3ch2n  -  Ifri2  +  nh2- 


hn  +  n2  +  h2 


’-h), 

where  ^CH.^  NH,  NH3>  N2>  and  H2  are  the  product  molecules. 


A  similar  scheme  may  be  written  for  MMH.  In  general,  we  wrice, 
for  the  hydrazine  family  of  fuels  (where  F  represents  the  fuel  molecule'' 


k 

hv  +  F-  A  +  NH2 
^1 

F  +  M - -  A  +  NH2  +  M 

^2a 

F  +  A  -£V  P1  -i  B 

W 

F  +  NH2  NH3  +  B 
^7 

E  +  NH2 - ►P1  +  N2  +  H2 


(v) 

(1) 

(2a) 

(2b) 

(7) 
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Tabulating  for  the  three  fuels  involved: 


where  P^  is  an  appropriate  product  molecule  and  A  and  B  appropriate  free 
radical  intermediates.  Where  more  complete  kinetic  representations  are 
required  for  more  complex  (analogous)  schemes  also  can  be  written 

for  UDMK  and  MMH. 

Clearly,  there  are  deficiencies  in  our  current  understanding  of 
the  kinetics  of  the  family  of  fuels.  The  calculational  techniques 

and  phenomenological  equations  discussed  in  the  text  are,  however,  quite 
general.  Thus,  the  general  features  deduced  in  the  body  of  this  work 
are  relatively  insensitive  to  modest  modifications  in  the  kinetics. 

This  is  in  fact  the  case  when  the  results  of  explosive  behaviors  of 
various  members  of  the  [^H^]  family  of  fuels  are  compared. 
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APPENDIX  II 

KINETICS  OF  DECOMPOSITION  OF  OZONE 

A.  OZONE  KINETICS 

A  rather  extensive  literature  exists  regarding  the  thermal  and 

A  2  6  7 

photochemical  propertied  of  czone.  The  simplified  kinetic  scheme. 

M+03**02  +  0+  M 
0  +  03  -  202 

represents  the  observed  data  quite  well,  for  dilute  ozone  in  the  absence 

4 

of  photochemical  sources.  For  concentrated  ozone  with  photochemical 
sources,  Semenov*^  has  argued  that  the  electronically  excited  states  of 

A  A 

0  and  02  (0  and  02,  respectively)  play  a  role  through  the  reactions 

* 

0^  +  hv  -*■  02  +  0 

°*  +  °3  *  °2  +  °2 
02  +  03  +  202  +  0 

Such  a  scheme  is  capable  of  providing  quantum  yields  greater  than  two. 

4  45 

Large  quantum  yields  (as  high  as  six)  have  been  reported  for  dry  ozone  * 

On  the  other  hand,  quantum  yields  of  less  than  two  are  also  reported, 
suggesting  that  any  usable  kinetic  scheme  must  also  include  such  steps 

38 

M  +  0*  ->•  M  +  0- 

t.  2 

M  +  02  +  0*  -*■  Oj  +  M 
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The  following  scheme  contains  chemical  kinetic  steps  that  are 
capable  of  representing  the  range  of  quantum  yields  observed,  as  well  as 
the  low  ozone  concentration  thermal  decomposition  data.  This  is  a  require¬ 
ment  of  an  acceptable  scheme. 


03  +  hv  — «L  02  +  0*  +  (Ey  -  ED  -  Ee) 

kl 

03  +  M  — W  02  +  0  +  M  +  h1 
k2 

03  +  0  — W  202  +  h2 

*  ^ 

o3  +  o  _ o2  +  o2  +  h3 

A  ^ 

03  +  02  —V  202  +  0  +  h4 

*  k5 

M  +  02  — ►  M  +  02  +  h5 
*  kfi 

M  +  02  +  0  03  +  M  +  h6 


(11-1) 
(II-2) 
(IX— 3) 
(II-4) 
(II-5) 
(11-6) 
(II— 7) 


In  the  following  section,  examination  of  this  scheme  is  made  with 
regard  to  quantum  yields  for  ozone. 


B.  QUANTUM  YIELDS  FOR  OZONE 

The  kinetic  scheme  given  by  equations  (1)  -  (7)  provides  the 
following  features: 

(a)  Thermal  initiation  as  well  as  photochemical  initiation. 

(b)  Collisional  deactivation  of  0_  (a  radiative  lifetime  greater 

-9  1 

than  10  sec  is  not  important  at  high  pressures) . 

(c)  The  possibility  that  quantum  yields  of  more  than  2  exist 
for  high  purity  ozone. 

(d)  The  possibility  that  quantum  yields  of  less  than  2  exist 
for  2ow  concentration  ozone  (in  02). 

(e)  The  p.  ssibility  that  quantum  yields  are  higher  at  high 
temperatures . 
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(f)  The  possibility  that  quantum  yields  are  pressure-dependent. 

(g)  The  possibility  that  quantum  yield«  are  not  independent  of 
the  radiative  flux  density  employed  in  their  experimental 
determination. 

/,  45  46 

Based  on  the  experimental  work  *  *  on  ozone  quantum  yields,  the 

aforementioned  features  appear  to  be  necessary  to  a  correct  representa¬ 
tion.  Examination  of  the  photochemical  steady-state  yields 
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The  quantum  yield  is  then  given  by: 
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Examination  of  Eq.  (11-12)  shows  that 


(II 


(a)  4>  decreases  with  increasing  flux  density  due  to  the  thermal 


initiation  step  and  the  fact  that  k  is  proportional  to  the 


0 
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2 

radiative  flux  density  (photons/cm '/sec) .  In  the  limit, 
where  k^  >>  k^[M],  <t  becomes  insensitive  to  k^.  In  the 
other  limit,  ky  <<  k^[M],  the  experimental  determination 
of  i>  is  of  questionable  value. 

(b)  For  high  purity  ozone,  [0^]  >>  [02]  and 


(c)  For  low  concentration  ozone,  [0^]  <<  [02]  and  <t  can  be  less 
than  2.  (For  very  low  concentration  ozone,  the  back  reac¬ 
tion  associated  with  Eq.  (2)  should  be  included.) 

(d)  For  high  purity  ozone  $  increases  with  increasing  temper¬ 
ature.  This  stems  from  the  fact  that  (k^/k,.)  increases 
with  increasing  temperature. 

(e)  For  high  purity  ozone,  the  temperature  dependence  of  k^,  as 
well  as  the  ratio  (k./k,.)  can  be  determined  from  the  quantum 
yield. 

In  order  to  employ  this  kinetic  scheme,  the  k's  must  be  known. 

For  a  given  case,  k  i3  determined  from  known  optical  extinction  coeffi- 

51  66  v 

cients  ’  and  the  radiative  flux  density  employed.  Values  for  k^,  k2, 
k^  and  k^  may  be  selected  from  those  reported  in  references  4,  55,  57, 

63,  and  64.  k,.  and  k^  may  be  deduced  from  quantum  yield  data.  In  parti¬ 

cular,  kj  may  be  deduced  directly  from  quantum  yield  data  for  high  con¬ 
centration,  high  purity  ozone  (once  k^  is  known) .  kft  is  then  deduced 
for  an  intermediate  case.  (e.g.,$  v  1). 


UNCLASSIFIED 


61 


UNCLASSIFIED 


APPENDIX  III 

CLASSICAL  ANALYSES  OF  THE  STABILITY  OF  THE 
PHOTOCHEMICAL  STEADY  STATE 

If  a  sample  of  an  unstable  compound  is  kept  at  a  bath  temperature 
below  that  required  for  explosion,  it  is  taken  to  be  in  a  quasi-steady 
state.  If  the  sample  is  then  subjected  to  a  thermally  and  photochemically 
significant  steady-state  radiative  field,  the  sample  wi1"  then  reach  a 
new  steady  state  (photochemical  steady  state)  or  explode.  In  general, 
the  steady-state  radiative  flux  density  necessary  to  cause  such  explo¬ 
sions  is  small  compared  to  the  pulsed  radiative  flux  densities  employed 
in  typical  flash  photolysis  experiments.  In  this  section,  sample  calcu¬ 
lations  are  made  using  thermal  theory  analyses  (generally  acceptable  at 
.high  pressures)  of  the  ordinary  and  photochemical  steady  states. 

A.  MIXED  HYDRAZINE  FAMILY  OF  REACTANTS 

For  a  cylindrical  sample  of  N^H^  gas  at  a  pressure  of  A6  atmospheres, 
a  bath  temperature  of  «  502. 5°K  implies  an  average  N^H^  temperature 
of  T^  *  524. 1°K.  These  are  the  critical  conditions  for  explosion  in  the 
absence  of  a  photodissociative  radiative  field.  If  the  sample  size  is 
made  smaller,  the  critical  temperatures  (T^  and  T^)  required  for  explo¬ 
sion  increase.  A  plot  of  the  steady  state  hydrazine  sample  temperatures 
is  given  in  Fig.  III-l,  for  three  different  sample  sizes.  The  steady 
state  temperature,  T^  ■  TCg»  corresponds  to  a  stable  sample  where 
(dTco/dT,  )  is  positive  and  finite.  Critical  explosion  conditions  corre- 
spond  to  (dTgg/dT^)  infinite  (T^  -  Tgg).  The  saraP~e  temperatures  for 
which  (dTgg/dT^)  is  less  than  zero  are  all  unstable  (Tgg  «  T2) .  The  data 
of  Fig.  III-l  are  for  zero  radiative  field. 

The  following  table  gives  the  variation  of  the  photochemical  steady 

state  properties  corresponding  to  the  same  sample  of  hydrazine  under  u.v. 

-2 

irradiation.  It  is  seen  that  very  low  radiative  flux  densities  (a.  10 
2 

watts/cm  )  are  adequate  to  assure  sizeable  [NHg]  concentrations. 
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d^  *  6  x  10  cm: 

Tcr1t  "  524-13  °K*  'fb’crU  '  502-5  'K 


dg  *  6  x  10"‘  cm: 

Tcrit  B  6S6‘05  °K>  (Tb)Cr1t  *  630‘31  °K 
dg  s  6  x  10"3  cm: 

TcrU  *  879‘06  °K’  (Tb}crit  *  833*5  °K 

600  700  800  900  1000 

ty  ck) 


Figure  III-l.  Steady-State  Reactant  (AVG)  Temperatures  for  100%  [N„H. ] 
at  P  ■  46  ATM  for  Various  (cylindrical)  diameters  and  4 
Bath  Temperatures. 
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Interestingly  enough,  ~  524. 1°K  is  again  the  critical  sample  tempera¬ 
ture  for  explosion.  For  a  given  bath  temperature  the  critical  [Nl^] 
concentration  is  specified  by  the  ci  .cical  radiative  flux  density. 

For  the  case  of  Monomethylhydrazine  gas  at  46  atmospheres,  thermal 
stability  calculations  yield  the  following  results: 


d 

cm 

T, 

be 

°K 

T 

12 

°K 

0.006 

666 

717 

0.060 

493 

520 

0.600 

394 

411 

For  the  case  of  unsymmetrical  dimethylhydrazine  gas  at  46  atmo- 

spheres,  thermal  stability  calculations  yield 

the  following  results: 

d 

T, 

T 

be 

12 

cm 

°K 

CK 

0.006 

583 

629 

0.060 

433 

457 

0.600 

346 

361 

B.  OZONE 

As  illustrated  for  the  mixed  hydrazine  family  of  reactants,  thermal 

stability  analyses  of  ozone, 

with  and  without 

radiative  fields,  have  been 

performed.  Results  for  gaseous  ozone  at  a  pressure  of  1  atmosphere  are 

given  in  the  following  table 

.  It  is  seen  tha 

t  rather  small  radiative 

flux  densities  have  large  effects  on  T^c  as  well  as  the  photochemical 


UNCLASSIFIED 


65 


7.5  *  :o° 


2.11  x  10 


1.11  x  10*1  j 
7.50  x  10"2  I 


3. *5  x  10 
1.97  x  10 


1.11  x  10 


2.1?  x  10 
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steady-state  concentrations  of  [0],  [0  ],  and  [02].  As  the  sample  size 
and  bath  temperature  are  decreased,  the  photochemical  contribution  to  the 
explosion  requirement  becomes  increasingly  more  important. 
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APPENDIX  IV 

EXPERIMENTAL  STUDIES  PERFORMED  ON  EXPLOSIONS  OF 
HYDRAZINE  AND  OF  OZONE 


The  purpose  of  the  experimental  program  was  to  demonstrate,  for  a 
few  selected  cases,  the  relative  roles  played  by  radiative  fields  in 
explosion.  In  support  of  the  DRC  program,  these  experiments  were  carried 
out  by  Rocketdyne  (A  Division  of  North  American  Aviation,  Inc.)  at  their 
Canoga  Park,  California  facility.  Before  presenting  the  text  of  the 
Rocketdyne  report  on  the  experimental  findings,  the  significance  of  these 
findings  will  be  noted. 

The  two  systems  studied  experimentally  are  hydrazine  and  ozone. 

.The  salient  conclusions  to  be  drawn  are: 

1.  It  is  far  easier  to  flash  initiate  an  explosion  in  ozone 
than  in  hydrazine.  This  experimental  finding  corresponds 
to  the  results  of  our  theoretical  investigation. 

2.  Hydrazine  becomes  thermally  unstable  just  below  300°C  and 
can  decompose  rapidly  even  where  a  sharp  temperature  or 
pressure  pulse  does  not  occur.  From  the  point  of  view  of 
explosion  theory,  such  thermal  instability  is  interpreted 
as  an  "explosion."  Time-dependent  examination  of  this 
thermal  instability  is  necessary  if  the  pressure-time 
and/or  tempsrature-time  history  of  the  sample  is  to  be 
described  Thus,  t^e  time-dependent  theoretical  studies 
performed  at  DRC  correspond  closely  to  the  experimental 
findings:  Hydrazine  may  become  thermally  unstable  (with 
or  without  radiative  sources)  and  can  decompose  to  pro¬ 
ducts  without  a  "sharp  explosion."  We  may  arbitrarily 
consider  a  "sharp  explosion"  to  be  characterized  by  a 
reactant  half-life  that  differs  by  no  nrre  than  some  5 
microseconds  from  the  reactant’s  three-quarter  life. 
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3.  Ozone  behaves,  explosively,  in  a  more  conventional  manner 
than  does  hydrazine.  The  experimental  results  indicate 
that  ozone  is  largely  unreacted  until  an  explosion  limit 
is  attained.  A  flash-initiated  explosion  occurs  rapidly 
(unlike  the  case  of  hydrazine) . 

The  unmodified  text  of  the  Rocketdyne  final  report  to  DRC  on  the 
experimental  work  follows: 

INTRODUCTION 

Experimental  studies  have  been  conducted  with  the  objective  of 
determining  conditions  for  photo- initiation  of  explosions  in  certain 
mono-propellant  compounds.  Since  the  Rocketdyne  study  represents  only 
a  portion  of  the  investigation  being  conducted  by  the  Defense  Research 
Corporation,  this  report  will  be  mainly  factual  with  the  theoretical 
interpretation  of  the  observed  results  appearing  later  in  the  formal 
report  for  the  entire  program. 

EXPERIMENTAL 

The  major  components  of  the  flash  apparatus  were  obtained  from 
Amglo  Corporation,  Chicago.  The  flash  lamps  used  had  a  linear  tubular 
quartz  envelope,  a  one  inch  arc  length,  and  were  filled  with  argon  and 
mercury  for  best  ultraviolet  output,  A  bank  of  ten  50-microfarad  con¬ 
densers  was  used  for  energy  storage  and  a  0-5000  volt  power  supply  was 
used  for  charging  the  condensers.  At  voltages  below  the  breakdown  vol¬ 
tage  of  the  lampe  (800  volts  for  lamps  mentioned  above),  a  trigger  supply 
which  provides  a  very  high  voltage  pulse  was  used  to  initiate  the  flash. 
At  voltages  above  the  self-flash  voltage  of  the  lamp,  a  special  relay 
was  used  to  switch  the  condensei  voltage  to  the  lamp,  which  then  flashed 
without  further  triggering.  For  most  of  the  work  reported  the  lamps  were 
used  at  their  maximum  rating  of  500  J  of  electrical  energy,  which  corre¬ 
sponds  to  500  microfarads  at  1400  volts. 
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An  inductance  of  600  u  H  was  used  in  series  with  the  lamp  to  reduce 
the  peak  current  and  thus  increase  lamp  lifetime.  From  information  pro¬ 
vided  by  the  supplier,  flash  duration  for  the  system  used  is  of  the  order 
of  400  microseconds  (to  1/2  maximum  intensity)  with  the  inductance  in 
the  circuit,  and  200  microseconds  without  the  added  inductance. 

Lamp  output  may  be  estimated  from  published  data.  Porter  reports 
20  0 

10  quanta  in  the  2000-4000A  interval  from  a  1000  J  flash  of  an  Xe  filled 
lamp  (Ref.  1).  Commercial  lamps  are  reported  to  have  efficiencies  of 
10-50  lm/W.s  (Ref.  2),  presumably  in  the  wavelength  region  of  approxi- 

O 

mately  3500-7500A.  Strong  (Ref.  3)  reports  that  their  Xe-filled  quartz 
lamps  produce  a  continuum  plus  Xe  lines  with  intensity  roughly  indepen- 

O 

dent  of  wavelength  down  to  3000A  (decreasing  somewhat  between  4600  and 

O  O 

5600A  to  a  minimum  at  5100A  of  about  60%  of  the  average  intensity) .  The 

O 

.wavelength  region  covered  was  2400-6500A.  Using  these  values,  it  is 
estimated  that  flash  lamps  used  in  the  present  study  produce  a  total  of 

O 

approximately  15  J  of  energy  as  light  in  each  1000A  intei  _i.  for  e  iOOO  J 
(electrical  energy  discharged)  flash  over  the  wavelengti  igion  of  2000 

O 

to  8000A. 

With  a  lamp  to  cell  distance  of  about  one  inch,  the  amount  of  light 
reaching  the  sample  is  estimated  to  be  about  1%  of  that  emitted  by  the 
flash,  because  of  geometric  factors,  scattering  by  the  quartz  sample 
tube,  etc.  This  estimate  assumes  a  6-mm  ID  sample  tube  one-inch  long 
with  a  reflector  behind  the  lamp. 

The  sample  tubes  employed  were  constructed  from  1  mm  ID  fused  silica 
capillary  and  3  mm  and  6  mm  ID  fused  silica  tubing.  By  sealing  water  in 
the  tubes  and  heating  slowly  until  the  tube  ruptured,  it  was  established 
that  the  1  mm  capillary  could  contain  pressure  in  excess  of  3000  psi;  the 
3  mm  tubing  ruptured  at  about  1000  psi  and  the  6  mm  about  400  psi.  The 
hydrazine  samples  were  sealed  in  sample  tubes  by  freezing  the  hydrazine, 
evacuating  the  system,  and  torching  the  sample  tube  shut.  The  ozone 
sample  tubes  were  filled  through  a  stopcock. 
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Heaters  were  constructed  of  aluminum  rod  drilled  to  contain  the 
sample  tube  and  a  thermocouple  and  wrapped  with  resistance  wire  and 
asbestos.  A  remote  control  device  was  constructed  to  insert  the  sample 
tube  into  the  heater  and  withdraw  It  to  a  position  next  to  the  flash 
lamp.  The  apparatus  was  surrounded  by  heavy  shielding  for  safety.  It 
was  found  that  on  insertion  of  the  sample  tube  into  the  heater,  the 
sample  reached  essentially  the  heater  temperature  in  two  minutes.  This 
method  of  operation  was  chosen  to  reduce  the  time  the  sample  was  heated 
prior  to  exposure  to  the  flash. 

Ozone  was  prepared  from  cylinder  oxygen  in  a  Pyrex  apparatus 
(Ref.  4)  consisting  of  two  tubes  sealed  with  an  annular  space  about  3  mm 
across,  where  the  oxygen  was  subjected  to  an  electric  discharge  from  a 
10,000  volt  ignition  transformer.  The  ozonu  °.r  was  operated  in  an  ice 
bath  at  0  C  to  improve  yields.  The  gas  mixture  was  analyzed  for  ozone 
by  bubbling  it  through  a  neutral  buffered  KI  solution  and  titrating  the 
I2  produced  (03  +  2H+  +  2I~  +  02  +  H20  +  I2>  with  thiosulfate  (Ref.  5). 

It  was  found  that  the  ozone  yield  at  atmospheric,  pressure  varied  from 
5.3  w/o  with  an  oxygen  flow  rate  of  130  cc/min  to  1  w/o  at  a  flow  rate 
of  830  cc/min.  The  ozone  was  handled  in  a  vacuum  system  with  stopcocks 
lubricated  with  Halocarbcn  grease  (Halocarbon  Products  Co.).  In  oper¬ 
ation,  a  170  cc  volume  was  flushed  out  and  filled  at  atmospheric  pressure 
with  approximately  57.  ozone  in  oxygen.  Then  liquid  nitrogen  was  placed 
on  the  attached  sample  tube  (volume  about  8  cc)  and  most  of  the  ozone 
and  oxygen  frozen  out.  With  the  sample  tube  in  liquid  nitrogen  the 
oxygen  was  pumped  off,  leaving  a  droplet  of  nearly  pure  ozone,  which 
was  easily  visible  to  the  eye.  The  stopcock  was  then  closed  and  the 
ozone  allowed  to  vaporize  in  the  8  cc  sample  tube.  For  higher  temper¬ 
ature  operation  the  sample  tube  was  heated  with  hot  air  blown  over  it 
and  the  temperature  was  read  from  a  thermocouple  in  a  similar  tube  next 
to  the  sample  tube.  Again  the  sample  apparatus  was  enclosed  in  shieldiag, 
and  heating,  operation  of  liquid  nitrogen  dewar,  etc.  were  carried  out 
remotely.  A  two-inch  thick  Lucite  block  allowed  observation. 
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RESULTS  AND  DISCUSSION 


HYDRAZINE 

When  the  hydrazine  sample  tubes  (3  mm  and  1  mm  ID)  initially  20% 
filled  with  liquid  hydrazine  were  heated  slowly,  they  ruptured  violently 
at  about  250  C.  It  is  believed  that  these  were  pressure  "explosions" 
caused  by  the  slow  thermal  decomposition  of  the  hydrazine.  Since  the 
critical  density  of  hydrazine  is  0.234  ya/cc,  the  tubes  could  not  have 
been  burst  by  hydrostatic  pressure.  The  vapor  pressure  of  hydrazine  is 
about  900  psi  at  290  C — below  the  measured  bursting  pressure  of  the 
sample  holder  tubing.  However,  the  calculated  final  pressure  of  non¬ 
condensible  gas  (N^  and  possibly  if  slow  decomposition  should  occur, 
4500  psi  at  300  C,  exceeds  the  burst  pressure  of  the  tubing. 

One.  mm  tubing  initially  about  5%  filled  with  liquid  hydrazine  did 
not  rupture  after  heating  at  300  C  for  several  hours.  At  the  conclusion 
of  the  heating,  liquid  was  visible  in  the  tube  (at  room  temperature)  in 
quantity  slightly  less  than  initially  present.  This  liquid  vaporized 
immediately  at  room  temperature  when  the  tube  was  opened,  leading  to 
the  conclusion  that  most  of  the  hydrazine  had  decomposed  to  form  ammonia, 
since  hydrazine  vaporizes  slowly  under  these  conditions.  The  calculated 
final  pressure  in  this  case,  1100  psi,  is  less  than  the  3000  psi  necessary 
to  burst  the  capillary  tubing. 

A  variety  of  samples  of  hydrazine  were  flashed,  but  since  no  explo¬ 
sions  occurred,  only  the  most  severe  conditions  will  be  mentioned.  In 
one  experiment  with  hydrazine,  an  apparent  explosion  did  occur  but  subse¬ 
quent  investigation  revealed  that  the  flash  lamp  probably  exploded  taking 
the  sample  cell  with  it.  The  inductance  was  placed  in  the  circuit  in  the 
remainder  of  the  experiments  to  greatly  reduce  the  probability  of  a  lamp 
failure  of  this  type. 

Liquid  hydrazine  was  flashed  at  room  temperature  in  a  3  mm  tube 
with  no  explosion.  6  mm  tubes  with  volume  about  1  cc  and  containing 
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0.05  gtn  hydrazine  were  placed  in  the  315°C  heater  for  2  minutes  and  then 
flashed  with  no  explosion.  These  tubes  later  ruptured  while  being  heated 
further  to  340-360  C.  The  results  indicace  than  an  explosion  did  not 
occur  and  would  have  been  detected  had  it  occurred. 

It  is  concluded  that  initiation  of  an  explosion  in  hydrazine  with 
light  is  relatively  difficult,  even  at  temperatures  at  which  thermal 
decomposition  is  occurring  rather  rapidly.  Because  of  this  thermal  de¬ 
composition,  attempts  at  higher  temperatures  are  not  feasible  unless  a 
very  fast  heating  technique  could  be  developed.  Since  it  appears  that 
hydrazine  must  be  heated  near  its  thermal  decomposition  temperature  before 
photoinitiated  explosions  can  be  observed,  it  is  recommended  that  in 
future  studies  the  experiment  be  designed  so  that  the  concentration  of 
the  hydrazine  can  be  monitored  continuously  up  to  the  time  the  flash  is 
.triggered.  The  inclusion  of  a  pressure  transducer  in  the  system  would 
also  ba  quite  useful. 

OZONE 

Calculations  based  on  the  oxygen  vapor  pressure  and  cell  volumes 
predict  an  ozone  pressure  of  600  mm  in  the  sample  tube  at  room  temper¬ 
ature.  This  was  doubled  in  some  cases  by  filling  the  measuring  volume 
twice.  No  difficulty  was  experienced  in  producing  or  handling  the  ozone. 
It  was  stable  for  24  hours  in  the  apparatus  at  room  temperature  as  evi¬ 
denced  by  visual  observation  of  readily  condensible  ozone  after  that 
time.  The  ozone  was  heated  to  70eC  for  a  period  of  about  a  minute  and 
was  again  observed  visually  on  freezing  down;  therefore,  thermal  decompo¬ 
sition  was  not  a  problem  under  the  conditions  employed  in  this  investi¬ 
gation.  Available  rate  constants  substantiate  this  conclusion  (Ref.  4). 

With  a  lamp  positioned  about  4.5  cm  from  a  6  mm  sample  tube  (no 
reflector) ,  ozone  was  still  observed  in  about  the  same  amount  on  freezing 
down  after  flashing  at  room  temperature  when  either  600  or  1200  mm  of 
ozone  were,  employed.  With  a  lamp  positioned  2  cm  from  the  sample  tube 
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and  with  an  aluminum  foil  reflector  1/2  inch  behind  '.ne  lamp,  no  ozone 
was  visible  on  freezing  down  after  one  flash  (500  J)  at  either  25  or 
70°C  and  with  600  or  1200  mm  pressure.  The  tube  did  not  rupture,  pre¬ 
sumably  because  of  the  relatively  low  pressure  of  gas  (even  if  the  condi¬ 
tions  of  the  decomposition  approached  adiabatic) .  Because  of  the  apparent 
disappearance  of  all  of  the  ozone,  it  may  be  assumed  that  a  very  rapid 
reaction  (explosion)  occurred.  McGrath  and  Korrish  (Ref.  6)  have  shown 
that  the  decay  of  ozone  after  a  flash  is  very  rapid,  and  the  presumption 
of  an  explosion  on  the  evidence  of  substantially  complete  destruction  of 
ozone  seems  to  be  valid. 

An  attempt  was  made  to  explode  condensed  ozone  by  flashing  imme¬ 
diately  upon  removal  of  the  liquid  nitrogen  bath,  before  the  ozone  evapor¬ 
ated.  There  was  no  apparent  effect  of  the  light. 
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This  concludes  the  unmodified  text  of  the  Rocketdyne  final 
report  to  DRC  on  the  experimental  work. 
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APPENDIX  V 


THE  COMPUTER  PROGRAM 


A  computer  program  (HYDRA)  was  written  in  FORTRAN  for  the  Control 
Data  fV>r  ation  3600  computer  to  perform  three  types  of  calculations  at 
the  option  of  the  user.  The  main  function  involved  the  integration  of 
Eqs.  (16)-(20)  of  Section  VI  to  describe  the  time  dependent  behavior  of 
the  system.  For  this  a  fourth-order  Rur.ge-Kutta  integration  with  vari¬ 
able  integration  interval  was  used.  The  remaining  two  options  allow  the 
U3er  to  iterate  on  the  bath  temperature  or  initial  radiation  flux  to 
determine  the  critical  value,  i.e.,  that  which  makes  R  and  L  tangent 
(see  Section  V) . 

The  input  consists  of  50  cards  for  the  first  case  after  which  one 
or  more  parameters  may  be  varied  for  subsequent  cases  by  adding  the  appro¬ 
priate  cards.  Specifically,  the  Input  is  as  follows: 

Card  1)  Fortran  format  3E20.10,  215 

DX:  Initial  integration  interval 

AE:  Absolute  allowable  error  for  integration  interval  control 

RE;  Relative  allowable  error  for  integration  interval  control 


IFVD: 


0  variable  interval  control 
1  fixed  interval  is  used  for  integration 


NMAX:  maximum  number  of  integration  steps  allowed 


Card  2)  Reactant  identifier  card 

Cols.  13-39  identifier  of  fuel  F  (e^g.,  ^H^) 

Cols.  23-29  identifier  of  reactant  A  (e.g.,  NH) 

Cols .  33-39  identifier  of  reactant  B  (e.g., 

Cols.  43-49  identifier  of  reactant  C  (e.g.,  N^H^) 
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Cards  3-49)  Parameter  List 

Cols.  1,  2  code  identifying  parameter 
Cols.  13-19  name  of  parameter 


Cols 

.  21- 

-40 

value 

of  parameter  in 

E20.ll 

format 

Cols 

.  55- 

-70 

units 

which  p, 

rrameter 

is  given  in 

According  to 

the  following  list 

• 

■ 

Code 

Parameter 

Code 

Parameter  Cods 

Parameter  Code 

Parameter 

01 

P 

0 

14 

h6 

27 

P 

o 

40 

\b 

02 

T 

'b 

15 

h7 

28 

P1 

41 

c4 

X 

03 

Fo 

16 

\ 

29 

SP 

*1 

42 

t 

0 

04 

p 

17 

\ 

30 

dS_> 

P1 

43 

T 

0 

05 

C 

P 

18 

\ 

31 

P2 

44 

F 

o 

06 

e 

19 

\ 

32 

S 

P2 

45 

A 

o 

07 

X 

20 

\ 

33 

dS 

P2 

46 

B 

o 

08 

d 

21 

\ 

34 

ATb 

47 

C 

o 

09 

ae 

22 

X 

35 

dTb 

99 

end 

10 

ah 

23 

X 

36 

A  P 

o 

11 

hI 

24 

X 

37 

dP 

o 

12 

h2 

25 

X 

38 

t 

max 

13 

h5 

26 

\t 

39 

Cdf 
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The  first  parameter  gives  the  pressure  which  is  not  used  in  the 
calculations  and  is  included  for  reference  only.  Parameters  2-5,  7,  8, 
10-15  are  as  in  Section  VI.  Parameter  6  and  26-33  are  used  to  give  the 
time  dependent  according  to  the  formula 

2  2 
k  -*  C,  (moles/watt-sec)  x  e(cm  /mole)  x  P  (watts/cm  ) 

v  v 

IP  watts/cm^  0  <_  t  <_  S_  -  dSp  /2 
o  if^ 

P.  watts/cm^  S  +  dSp  /2  £  t  <_  Sp  -  dSp  /2 
1  P1  *2  *2 

P„  watts/cm^  Sp  +  dS  !2  <_t 
/  r  2  if  2 

In  the  transition  regions,  Sp  -  dSp  /2  _<  t  <_  S  +  dSp  /2  ; 

*1  1  1  *1 

.Sp  -  dSp  /2  <_  t  £  Sp  4  dSp  /2,  the  value  of  P  changes  smoothly  from 
P2  *2  *2  *2 

one  value  to  the  next  using  a  fifth  order  polynomial  having  zero  first 

and  second  derivatives  at  the  end  points.  The  value  of  C,  (parameter 

vT 

26)  is  used  to  compute  k  in  dT/dt  while  C,  (parameter  40)  is  used  to 

V  K  « 

vB 

compute  kv  as  it  enters  into  dB/dt.  This  allow3  the  user  to  distribute 
the  effect  of  the  radiation  into  thermal  energy  or  dissociation  in  any 
ratio  he  chooses. 


Parameters  16-25  are  used  to  compute  the  rate  constants  according 
to  the  formula 


k 


i 


<\i  ■  VR> 


Parameters  34  and  35  are  used  in  the  iterations  of  to  determine 
the  critical  value.  One  must  input  a  value  of  T^  as  a  negative  number, 
for  which  the  R  and  L  curves  intersect.  Then  the  program  starts  increas 
ing  Tfe  by  the  amount  AT^  until  this  condition  is  violated.  It  then  pro¬ 
ceeds  to  converge  on  the  critical  value  until  consecutive  values  of  T^ 
are  within  the  range  i*.  . 
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Parameters  AP  and  dP  perform  the  same  function  for  the  k  iter- 
o  o  vo 

ations  and  the  same  rules  apply. 

In  some  cas^.s,  the  user  may  want  the  time  dependent  behavior  of  the 
system  under  conditions  where  the  fuel  concentration  remains  constant. 
Parameter  39  performs  this  function  as  it  multiplies  the  expression  for 
dF/dt  in  the  program.  Thus  it  should  have  the  value  1  for  normal  oper¬ 
ation,  0  when  the  fuel  concentration  is  to  remain  .opstant. 

Parameter  38  acts  as  a  control  on  computer  time.  It  represents 
the  maximum  value  of  t  to  which  integration  can  proceed. 

Normally,  the  initial  conditions  on  the  dependent  variables  (with 
the  exception  of  F  )  are  taken  a.s  the  steady  state  values.  Parameter  41 
overrides  this  if  it  is  input  as  any  non-zero  quantify.  In  that  case 
the  initial  conditions  are  taken  from  parameters  42-47. 

The  last  card  merely  indicates  to  the  program  that  it  has  completed 
reading  data  for  a  particular  case. 

The  user  has  the  option  of  using  his  own  identification  for  the 
above  parameter?  (1-47);  however?  some  identification  should  be  entered 
as  the  program  writes  this  same  information  as  output  to  identify  the  '• 
case. 


If  one  wishes  to  run  additional  cases,  only  those  parameters  which 
change  need  additional  input  cards  followed  by  a  "99"  card. 

.-ddltional  remarks  on  the  data  contained  ir.  the  first  card 
should  be  made.  The  integration  subroutine  used  in  the  program  controls 
the  integration  interval  size  by  comparison  of  Runge-Kutta  integration 
over  two  half-intervals  with  a  Simpson’s  rule  sum  of  the  resulting  values. 
Thus  if  AICD  is  the  incremental  change  in  the  dependent  variable  according 

bK 


UNCLASSIFIED 


to  Simpson' 8  rule  and  AI^,  is  the  same  as  determined  by  Runge-Kutta,  the 
criterion  for  increasing  or  decreasing  the  interval  is  determined  by  the 
term 

airk  ,  -aisr 

A  +  RI 

where  I  is  the  current  value  of  the  dependent  variable  and  A  and  R  are 
the  absolute  and  relative  allowable  errors  given  on  the  first  card. 

The  program  listing  follows  with  the  exception  of  the  RKS3  numerical 
integration  subroutine.  This  is  a  standard  subroutine  available,  for 
example,  through  the  IBM  Share  Library. 


»;  7.'”'-'  y&r  -v  r  i.?  ■  ’•  yrr" 

mr--- ' ' 

•/!>  t- 


j  -»  ,  *  \*<*  ■u‘«V,<»T' 
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PROGRAM  HYDRA 

COMMON  Z# LABEL# DATA# NAME# ZS# US# TS#DR$#XLS#Ci#C2,C3#C4, 
t  CS#C6#C7, 

X  MS# C8 »  TE S  Tl#  TEST2#  D J  V #  ULi DUL#  UST ART # UCR1#UCR2# TCR1# 
t  TCR2#RPR#XLPR# 

X  F,EI,T1#T2 

COMMON/ 1/  I XI#  IK2# |K5#1K6# IK/# | PUMP , NCOUNT . NHAX 
COMMUN/OE/ 

♦U.DU#S#DS,30RK#AE»RE#T3#T4#T5#T6#T7#T8»T9«VP  #SS*SCH1» 
X  SCH2#SCK3,5CH4,  AjL*  A2#P#T #  DEI# DK1 # DE2# DK2 # DE5# DK9 # D£6# 
t  DK6#DE7*DK7, 

X  TPU»C9#C10,ST0P#FSTOP 
DIMENSION 

t  Z<10(n*LABEL<loa>*DATA{8)»NAME(4),ZS(2)#TS(2)#US(2)« 
X  RS(2),DRS(2)#XLS<2), UNITS! 100,2) 

DIMENSION  U(5)iDJ(5)#AE(&)#RE(5)#woRK(50I#VP(51,13) 
EXTERNAL  PHD#PHC 

2NTEQRATION 

♦  PARAMETERS  ON  THE 

FIRST  CARD 


1C0U  READ  ( 60 « 1010  >  DS*AO,RO, JFVD#NKAX 
1010  FORMAT! 3E20«li*2l5) 

DO  1011  1«1,10Q 
10U  Z(I)  *  0,0 


READ  REACTANT  NAMES 

AND  INPUT 

PARAMETERS 


i02u  in  •  x 

1029  READ ( 60 # 1030 )  INPiDATA 

1039  F0RMATU2#8X#7A8*A6) 

1040  IF  ! EOF  #  6 U  }  1041,101)0 

1041  CALL  HXU 

1050  IF  1 |NP,ESS,99)  BO  TO  1121 
1060  IF(JNP,NE,0)  GO  TO  U00 
1070  D6C0UE<64, 1080, DATA)  NAME 
If  80  F0RMAT!4(2X#R8>) 

1090  GO  TO  1020 

MOO  DECODE <64,1110#DATA) 

♦  LABELUNP)#ZUNP)#(UNITSCINP#J)9  J»l#2) 

1U0  FORMA7'2X#A8#10r  E20 , 11  #  4X#  2AB ) 
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1120  GO  TO  1020 

C  WRJTS  PARAMETER  LIST 

C 

1121  lF(Z(a;.QT.OtO)  GO  TO  1122 
Z(2)  s  -Z<2> 

nr  *  2 

UETB  «  Z(34) 

CONTB  *  2(35) 

DIVT  «  1,0 
GO  TO  1130 

1122  JF(Z(27).GE,Q.0>  GO  TO  1130 
2(27)  «  •Z<27> 

1 1  T  *  3 
DfeP  •  HU) 

CONP  «  Z (37) 

UlVP  *  1,0 

1130  WRITE(6i»1140)(NAMECi)| I«1,4>,(J»LA3EU(J)#Z(J) 

f  .  ,(UN!TS(J,K),X«1»2> 

X  » J«l»47) 

1140  F ORMAT ( 1H1 , 45H  MIXED  HYDRAZINE  FUELS  •  EXPLOSIONS 

f  ,//// 
X  20 X  THE  REACTANTS  ARE  //  6H  F  =  A8*6H  A  *  A8.6K 
r  B  *  A8»6H  C 

X  *  hH,//  20M  THE  parameter  list  *// 

f  <2X, J2,2X,A3#5Hf , , . . , E20 , 11, 4X 

X  ,2A8>) 

C  COMPUTE  CONSTANTS 

C 

115U  Cl  •  Z(9)/(2,0*Z(10))  ♦  1,0 

1160  C2  *  Z(10)*Z(3)*(Z(H)**2)/(8,0*Z(7) ) 

1170  IF ( Z ( 16 ) , EO , 0 , 0 )  GO  TO  1210 
1180  C3  4  LOGF ( Z ( 16 ) ) 

1190  1 K 1  *  1 
1200  GO  TO  1220 
121U  I K 1  *  2 

1220  IF ( Z ( 18) |EQ • 9 , 0  }  GO  TO  1260 
1230  C4  •  LOGF (Z(18) ) 

1240  IK2  «  1 
1250  GO  TO  1270 
1260  IK2  «  2 

.1270  .(.F  C Z € 20*1 , EO.O  « 0  )  GO  TO  1319 
128U  C5*L0GF (Z(20 ) ) 
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1291)  IKS  *  1 
1300  GO  TO  1320 
1310  t«5  a  2 

1320  IF ( Z ( ?2  5  ,  EO  .  0 , 0  5  GO  TO  1370 
1340  C6  *  L0GFU<22)) 

1350  IK6  *  1  * 

1360  GO  TO  1360 
1370  I K6  a  2 

1380  IF ( Z ( 24 )  | EO • 0 « 0 )  GO  TO  1420 
1390  C7  a  LOOK (Z (24)) 

1400  IK7  «  1 
1410  GO  TO  1430 
1420  l K7  3  2 

1430  08  a  l6,Q*Z<7>/</.<8)**2> 

1431  STOP  *  10 1 Q*Z ( 2 ) 

FSTOP  3  U|1*Z<3) 

1432  C9  *  Z ( 6) 

CIO  3  C8/ ( Z ( 4 } *Z ( 5 ) ) 


C 

c  TEST  FOR  STABILITY 

C 


1440  TEST1  3  UR( Z ( 2 ) )  »C8 
1450  iF(TfeSTlilT.Q.O)  UU  TO  1510 
1460  TEST2  3  H(Z(2)) 

1470  lF(TfeST2.l,T. 0.0)  GO  TO  1500 

1480  GO  TO  (1483,1481#1481)  I  XT  - 

1481  WRITt(61il482)  Z(2)»Z(27) 

<182  FORMAT  ( //  6H  Tt?  »  E20  ,  It#  I0X.  6H  PO  »  E20,ll  ) 

1483  HH I Tfe( 61*1490 )  T6ST2. TEST1 

1490  FORMAT ( ///64H  IMPOSSIBLE  TO  MEET  STABILITY  CRITERIA 

t  AT  TB  R-L  HA 

XS  THE  VALUE  E20.U.//  26H  THE  DERIVATIVE  OF  R*L 

♦  IS  E20-.ll) 


1500  GO  TG(1020#2490<260U)  II T 

1510  DIV  t  jl.Q 

1520  UL  *  5.0 

1530  UUL  a  5,0 

1540  TL  «  UL  *  Z<2) 

1550  TESTl  *  UR(TL)  -08 
1560  IF  (  TfeST.l  >1570  #  1'640  ♦  1590 
1570  UL  «  UL  ♦DUL/DIV 
1580  GO  TO  1610 


-"V 
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1590  0 1 V  *  2,0 

1600  01  *  UL-DUL/DIV 

1610  lF(UUL'l£,l,OE-Oii)  GO  TO  1640 

1620  DUL  *  DUL/DIV 

1630  GO  TO  1540 

1640  USTART  a  UL 

1650  T6ST2  a  R(  1 L)*C8*UL 

1660  IF(TEST2)1/20i1700«1670 

1670  GO  TO  (1673, 1671, 1671)  JIT 

1671  WRJTE(61il402)  Z(2)#Z(27) 

1673  WR 1  Tfe ( 61# 168Q )  UL, TL# TEST2# TEST1 

1680  FORMAT ( ///  71H  IMPOSSIBLE  TO  MEET  STABILITY 

♦  CRITERION,  MINIMUM  QF 
X  R-L  OCCURS  AT  U  *  E20.U.10X#  5H  T  *  E2QtU»//  30H 

♦  WHERE  JT  HAS 
X  THE  VALUE  H*l  *  £2Q ,11, 10X# 10H  D(R-L>  •  E20.ll) 

1690  GO  TO  (1020,2490,2600)  JIT 

1/00  TCR1  a  TCR2  *  TL 

1701  UCRl  a  UCR2  «  UL 

1710  GO  TO  2020 

1720  0 1 V  a  S.Q 

1730  DUL  «  0,1*UL 

1740  GO  TO  1780 

1750  TL  «  UL*Z(2) 

1760  IEST2  a  H ( I L )  -C8*UL 

1761  1F(UEST2,LE,0,0),AND,  ( ABSF  (UL>  ,LE  •  l .  06-08  )  )G0  TO  1850 

1770  IF ( TfeST2 ) 1780 • 1850,1800 

1771  UCRl  a  UL 

1780  OL  «  MAX1F(0.0, UL-DUL/DIV) 

1790  GO  TO  1820 

1800  OJV  a  2,0 

1801  JCN1  a  UL 

1810  UL  «  UL+DUL/O JV 

1820  IF ( DUL ,LE. 1, 0E-08)  GO  TO  1860 

1830  DUL  *  DUL/DIV 

1840  GO  TO  1750 

1850  UCRl  a  UL 

1860  TCRl  a  UCRl  *Z(2) 

1870  DUL  «  0,1*(USTART) 

I860  UL  *  USTART  ♦  DUL 
1890  DlV  a  1,0 
1900  fl  *  UL  ♦  2(2) 


pwyjtmmBF^.M.wy.^w  wwww^Bjwppww 
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1910  1EST2  s  R(1L)  -  C8*OL 

1911  IF((TEST2,LE, 0,0), AND, <TL,QT, STOP))  QO  TO  2000 
1920  IFITEST2J193Q, 2000,1960 

1930  UCR2  *  Ul 

1940  Ul  *  UL  +  UUL/DIV 

1950  QO  TU  1990 

I960  DIV  *  2,8 

1970  UCR2  *  Ul 

1980  Ul  «  Ul-DUL/DIV 

1990  1F(DUL,LE,1,OE-08>  00  TO  2010 

1991  DUL  *  UUL/DIV 

1992  QO  TO  1900 

2000  UCH2  *  UL 

2010  TCR2  *  UCR2  *Z<2< 

2020  UL  *  UCR1 

2030  DUL  *  ( UCR2-UCRD/4Q,  0 

2040  IF ( DUL , EU , 0 • 0 )  GD  TO  2132 
QO  TU  (2041,2043,2043)  IIT 

2041  wH I T£( 61, 2042 ) 

2042  FORMAT ( 1H1, 53H  THE  FUNCTIONS  R  AND  L  BETWEEN  POINTS 

♦  OF  INTERSECT!1! 
XN  //  14X,1HU,24X.1HT.24X»1HR.24X,1HL,20X,4H  R*L  //) 

60  TO  2050 

2043  WR I TE ( 61 , 2044 )  Z(2)*Z(27) 

2044  FORMAT { lMl, 67H  THE  FUNCTIONS  R  AND  «.  BETWEEN  THE 

♦  ®OINTS  OF  INTERS6 
XCTION  FOR  re  *  E20 ,11,10X#5HP0  *  £20,11,// 

♦  14X,1HU,24X,1HT,24X,1HR 
X,24X,1HL,20X,4H  R-L  //) 

2050  DO  2130  1=1,41 

2060  TL  «  UL  *Z(2) 

2070  RPH  *  H  C  r L  ) 

2080  XLPR  *  C8*UL 

2090  TEST2  *  RPH  -XLPR 

2100  wR 1 Tfe( 61, 2110 >  Ul,TL,RPR,XLPR,TEST2 

2110  FORMAT (  5( 4X,E20 , 11 ) ) 

2120  UL  «  UL  ♦DUL 

2130  CONTINUE 

2131  QO  TO  2140 

2132  WR1TE(61,2133)  UCRA 

2133  FOHMaT(///  32H  R.  AND  l,  ARE  EQUAL  ONLY  AT  U  *  E20.il) 

C 
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C  FIND  TS  AND  TSP 

C 

2140  2SU)  *  0,0 
2150  ZS<2)  a  Z(26) 

2160  UO  2380  K*li  2 
2170  DIV  a  1,0 

2180  ul,  *  0,0 

2190  UUL  •  5,0 
2200  Z  <  26 )  a  ZSUG 
2210  TL  «  UL  ♦  2(2) 

2220  F  a  C8*UL  -  R(TU 
2230  IF(F)2240, 2330*2270 
2240  US ( K )  a  UL 
2250  UL  «  UL  *  DUL/iiiV 
2260  UO  TU  2300 
2270  DIV  •  2,0 
2380  USIK)  a  UL 
,2290  UL  a  UL-DUL/D1V 
2300  IF  <  DUL  , l.fe.  1 . 0E*08)  GO  TO  2330 
2410  DUL  *  DUL/OIV 
2320  Uti  TU  2210 
2330  US(K)  a  UL 
23«0  TSIK)  a  US(K)  *Z(2) 

2350  RSlK)  *  H(rS(K)) 

2360  DRS(K)  a  DR( TS(K ) ) 

2370  XLS(K)  a  C8*UL 
2380  CONTINUE 
2390  WRlTE(61f 2400) 

t  US<l)iUS(2),TS<l)#TS<2)#HS<i)#RS<2>»XLS<l)iXL$(2)i 

X  DRS ( 1 ) , DRS(2 ) f C8f  C8 

2400  FOHKAT(///  4X.51H  THERMAL  STEADY  STATE  // 

♦  PHOTOCHEMICAL  STEADY  ST 

XATE  //  1X|2H  U»2X#E20,11#10X#E20,11»/  3H 
t  T,2X#520 ,11»10X|62Q ,11*/ 

X  3H  R, 2X*E20 , ill 10X, E20 i Ilf /  3H 

t  L ^  2 X i E2 0 i 11 1 1 0  X i E2 0  « 11 1 /  3H  DR# 

X  2X,£20,llil0X,E20,ll«/  3H  DL 1 2X  *  E3  0 , 1 1 » 1 0  X  *  £2  Q  j 11 ) 
2480  GO  TO  <2701 f 2520 f 2640 )  MT 
2490  DIVT  a  2,0 

2500  THH  a  Z<2) 

2501  Z<2)  a  Z(2)  -DET8/UIVT 
2510  GO  TO  2540 


cs  n 
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2520  TBL  *  Z(2) 

2530  HI)  *  Z(2>  *  DETB/DJ V? 

2540  iF ( DETB , LE . CQNTB )  GO  TO  2570 
2550  UETB  *  DETB/DJVT 
2560  GO  TO  1440 
2570  WRITE(61«258Q)  TBUTBH 

258U  FORMAT ( ///  37H  THE  LIMITING  VALUE  OF  fB  |S  BETWEEN 

♦  E20  till  6H  AND 
X  £20  1 11  ) 

2590  GO  TO  1020 
2600  UIVP  «  2,0 
2610  pH  *  Z<27> 

2620  Z<27)  «  2(27}  -  DfcP/DjyP 
2630  GO  TO  2660 
2640  PL  *  Z ( 27 ) 

2650  Z(27)  «  Z <27 }  ♦  UEP/DIVP 
2660  1F(DEP,LE,C0NP)  GO  TO  2690 
2670  OEP  *  DEP/DIVP 
2660  GO  TO  1440 
2690  WHITE(61i2700)  PL*PH 

2700  FOHHA7 ( ///  37H  THE  LIMITING  VALUE  OF  PQ  (S  BETHiEN 

♦  620  til#  6H  AND 
X  E20.ll  J 

GO  TU  1020 

SET  UP  INITIAL 

CONDITIONS  AND 
C  INTEGRATION 

♦  PARAMETERS 

2701  (F ( Z (41 ) ,NE«  9 « 0)  GO  TO  2741 
U<1)  *  TSM) 

U(2}  «  Z ( 3 ) 

NCOUNT  «  0 
U(3)  *  0,0 

IF((1K1,ED.2>,0R.(JK2,E0(2M  GO  TO  2710 
El  «  03*04  >(ZU7}*Z(19)>/(2,0*U(1)) 

(F (EltLEt*700,0)  GO  TO  2710 
T1  *  EXPF(El) 

GO  TO  2720 
2710  T1  *  0,0 

2720  lF(UK5,fc0.2},0R,(Iv<i7<EQtJZU  GO  TO  2730 
E2  «.C5-C7  *(Z(2iI*Z(25)>/(2,0*U(l)} 
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lF(E2,LE,-;00,0)  go  TO  273Q 
t?  «  expfte2) 

GO  TO  2740 
2730  T2  «  o,0 

2740  U  ( 4 )  m  Z ( 3 ) • ( Tl  ♦TZ > 

U(5)  «  o,o 

ilU)  «  US(1) 

S  «  0,0 
GO  TO  2743 

2741  S  «  2(42) 

UO  2742  J«l,5 

2742  U(J>  *  Z(42*J) 

2743  NPR  i  jj 
IBKP  a  l 

• NTHY  c  x 
ME  *  5 

00  2750  I*1#NE 
AEU)  «  AO 
27 5U  RE (  I)  3  ro 

If(Z(30) ,EQ,0,0)  GO  TO  2760 
bCHl  3  Z(29)  •Z(30)/2,0 
SCK2  .  Z(29)  ♦  ZcSo,/2?o 

SCH4  3  Z ( 32)  ♦  Z ( 33) 72 , 0  ' 

B0  T°  2,60 

2/6V;1U  HKSJIP““'PHC'UlOu'*e>HE.NORK.S,OS,Ne,irv0.1BKP, 

1820  "'"•IS*** 

*'n  "u‘mSRH17H  ERR0R  RETURN  0f  PROH  RKSJ  , 

GO  TO  1020 
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PRINTOUT  OF  A  TYPICAL  COMUPTER  RUN 
OF  THE  EXPLOSIVE  TRAJECTORY  OF  AN  IRRADIATED  UDMH  SAMPLE 
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-  ■  *  >  yjv  ,<rr  _ 1 


i.'iy'  IHI 
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•v  ft«N  334Mftft  ft  3  >  9  3  ft  ^  N  4)  9  ft  ®  4  4ft 

t  4O«4'»3*^4)Nft43B-J»N^t'3OOftf43-)4t!')>'ftft4H40X>9'J4£3tl4  ^  O -)  ft 
A»oiftftiftrwiywfliMOHO0«4(HO4fti3o©®noiV94nc»  ftcoftftoisiftM^oo 
n  *40NftW4»O00ftnftft"»3ft#*N''C»<H>t*#ft»4K0N'4N*N4*.<ON1ft44  9ftN 
>f>4»4tWN34NO4#9  3»Mft*H*)4»»^N»)4ftftftftft4«).\H>*«'0f4*ftM>ftH 
«  »93ft  «93ft»3ftlfl«N90H34  ftN  9>3»4«J4  ft«**.0»3H!W^4i»ft0>0>>  O  -4  *0  ft  X 
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